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O 

, The strong-coupling perturbation theory (SCPT) for correlated electron systems is extended 

to the case of full Coulomb interaction. The Coulomb mechanism of the orbital polarization is 
discussed and attention is paid to the importance of spectral weight transfer between the localized 
and delocalized subsystems of electrons. A one-to-one correspondence between subsets of Feynman 
^—K ' graphs of SCPT (which we name the approximation of renormalized Fermions (ARF)) and weak- 

ly^ , coupling perturbation theory (WCPT) is established. The comparison of the Galitskii and Migdal 

expression for the total energy and the Sham equation (which connects the self-energy and the 
exchange-correlation potential in density functional theory) for WCPT with the ones for the systems 
with strong electron correlations is used for the formulation of a simple theory for extending the 
D ■ local density approximation to density functional theory (DFT) to include explicitly correlations. 

The extension requires the inclusion of the many-electron spectral weights in the definition of the 
charge density and the renormalization of the mixing and hopping matrix elements caused by many 
body effects. 
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O I- INTRODUCTION 

The standard model of the lanthanides successfully explains most of magnetic properties and, within ab initio 
calculations, also the equilibrium lattice parameter, bulk moduli and cohesive energy; actually, the picture of localized 
/-electrons works much better for compounds with rare earths than the picture of fully delocalized /-electrons. 
Nevertheless, it is not sufficient for an explanation of all the complex physical and chemical properties. A recent 

■ example is given in the neutron scattering experiment on Pr metafi. The mechanism which can provide the /- 
localization, has been recognized long ago as the intra shell Coulomb interaction, which forms large energy gaps 

I . between .occupied and unoccupied d- or /-orbitals. It has been explained within Hartree-Fock approximation by 
BrandowB i±i connection with physics of Mott insulators and explicitly exploited in such phenomenological theories, 
as LDA+UlI, self-interaction corrected LDA DFTcrtl It is clear that this picture is a bit oversimplified, since in 

■ many cases it is difficult to describe a local subsystem of materials within a single-electron approach. LDA-DFT is 
not a single-electron theory since it goes beyond the Hartree-Fock approximation at least up to the random phase 
approximation (RPA). However, it is also well-known that without formation of an energy gap with the help of SIC 
or LDA+U between /(ei)-orbitals it is difficult to obtain electron localization in practical implementations of DFT. 

However, in total energy calculations within the local density approximation to density functional theory, it is 
possible to treat the /-electrons either as core (or valence) electrons, obtaining localized states. Irrespective of if the 
Q . /-states are localized or delocalized, there is a problem of taking into account strong correlations between /-electrons 
of non-filled shells, since the expression for the exchange-correlation potential usually is based on the expression 
. i-H , derived from the theory for homogeneous electron gas. The treatment of /-electrons as core states also requires a 
prescription of how many of them should be fotped to be localized in the core. Even though the total energy may may 
be minimized with respect to the /-occupations this approach must rely on experimental data of the /-shell and it has 
problems of including fluctuations or non-integral occupations of the /-states. A theoretical method that attempts to 
describe all the intricate properties of /- (and d-) electron materials, such as heavy Fermion behavior, Kondo effects, 
Mott insulators (with a transition from a metal to insulator), unconventional super-conductivity, hybridization gaps 
etc., must be build on a many body framework, and it seems that a fruitful way to do this is to join many body 
corrections and electronic structure methods. The main motivation for this combination is that the many body theory 
described accurately the correlated /- (sometimes d-) states, whereas electronic structure methods treat the diffuse, 
more extended states well. 

Here we suggest to look at the problem using a theory of strongly correlated electrons (SCE). We would like 
to combine field-theoretical methods with DFT. For the electrons which are localized the approach starting from 
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the atomic limit seems to be appropriate. The derivation of the connection, between DFT and many-body theory 
for the total energy given by ShamO uses the expression for the total energya which contains the sum of all skeleton 
graphs. The theory of SCE from the atomic limit is based mainly on the Pauli-principle (e.g. for destruction operators 
° 2 faa = 0)- Since some infinite sequence of graphs of perturbation theory from the weak-coupling side (say, for the 
self-energy for the Fermion Green function (GF)) corresponds to the expression where the relation like nf a(T = fif aa 
has been used within the strong-coupling perturbation theory, it is not easy to find this sequence and provide a 
one-to-one correspondence between these perturbation theories. It is possible, however, to express the Fermion GFs 
in terms of GFs for many-electron operators. Here we will exploit this connection. 

The description within the framework of the LDA to DFT, when the /-electrons are treated as core electrons with 
only a certain number of /-orbitals occupied may be provided within the many-body approach in the following way. 
Let us consider an ion which has n /-electrons in the ground state number. Then, only the transitions r„ — * T n ±\ 
will be involved in the formation of the spectrum of the single- electron excitations while all other transitions, like 
r n — > r n ±2, r„±3, involving larger number of electrons will be strongly suppressed by a large energy separation between 
these states. If we, in such theory, take the limit when the energy of the atomic-like transition A2 = E^ +1 ' — E^f 1 
between any n + 1 and n-electron states T and T' of the /-ion is much higher than the Fermi energy, e p , the number 
of /-electrons in the ion will be fixed. Indeed, in this limit this upper "single-electron" level A2 is empty. In the rare 
earth elements the populated part of the /-spectral density corresponding to the transitions Ai = E^) — £^"7 is 
much below ep. It can also be (much) below the bottom of the conduction-electron bands. This mechanism leads 
to a self-consistent switching off of the mixing interaction (see figures 4,5,6 for LDA+Hubbard-Kanamori corrections 
in Ref.cl) and removes the overlap between these core-like levels and the conduction electrons. Thus, this physical 
picture exactly corresponds to a type of ab initio calculation with a fixed number of localized electrons. The photo- 
electron spectroscopy experiments, show that even in the rare earth elements, for which this picture seems to be 
most appropriate, the level A2 is sometimes only slightly above Fermi energy. Therefore, due to mixing interaction it 
should contribute to the cohesive energyll3. The total energy of the pistem and Fermion GFs can be calculated in a 
simple approximation by making use of the theory developed in RefEH. This allows to construct the Sham equation 
for the exchange-correlation potential and to derive corrections to the standard LDA expression. 

The paper is organized as follows. In Section ||, in order to provide an opportunity to directly compare the 
contributions from Coulomb interaction to the equations of motion for GFs in WCPT a nd SC PT, the equations 
of moti on a re written exactly via functional derivatives of the GFs in WCPT (see Section II A ) and in SCPT (see 
Section II B ) within a real-space non-orthogonal representation. Here we introduce the approximation of renormalized 
Fermions (ARF), which allows us to establish a one-to-one correspondence between certain subsequences of graphs in 
WCPT and SCPT. This fact is used-further for the formulation of an recipe for a correlation corrected form of LDA 
on the basi s of the analysis by ShamQ for the exchange-correlation potential for a normal and a correlated system. In 
Section HI we discuss the standard model for the lanthanides and compare the self-energies within WCPT and SCPT 
in the ARF. Since the graphs which are not taken into account by this approximation can be written explicitly, this 
allows us to establish the domain of validity of the standard approach to the case when the system has (quasi)localized 
electrons, and, also, suggests a possible recipe to overcome the difficulties existing in the LDA in DFT. Section IV 
contains the conclusions and a discussion. The contributions to the equations of motion for the GFs are rewritten via 
functional derivatives of GFs. This allows to construct a regular SCPT via an iteration procedure for the equations 
for the GFs. In Appendix |b] we discuss the transfer of spectral weight between the low and high energy regions and 
show the role played by spectral weights in the formation of the orbital polarization in the simple example of 3-orbital 
atoms with 2 strongly interacting electrons. 



II. COULOMB INTERACTION AND GREEN FUNCTIONS 



A. WCPT for an Non-Orthogonal Set 

Perturbation theory from the atomic limit should, on one hand, be constructed in real space in order to be able to 
treat the single-site Coulomb interactions better than within the LDA. The wave functions centered on different atoms 
are not supposed to be orthogonal to each other. On the other hand, we have to compare the perturbation series for 
the GFs within weak coupling and strong coupling. For this reason we have to understand how the equations reflect 
the non-orthogonality of the basis set. Here, we shortly repeat the well-known Kadanoff and Baym derivation for the 
exact expression for the self-energy and obtain the first corrections which we will compare with the ones obtained 
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from the SCPT approach later. First we have to introduce the S'-matrix with an external field E/23(i) a 2 a 3^!]- Then, 
for the Hamiltonian 

W = h 23 a l a 3 + ^ w 234544 a 4 a 5: i 1 ) 

2 

where h^ 3 — (2| 9^7 + w(r)|3) and D2345 is a matrix element of the Coulomb interaction, the equation of motion for the 

GF Fn'(tt') = l/i(7~ai(t)a\,(t!)) has the form: (the letter F will be used for Fermion GFs below keeping G for the 
GFs involving Hubbard operators) 

{S 13 *d t - 0^[h° 23 (t) - <23]}F 31 ,(t,t'\U) = iS(t - t')0-?, 

+°r2%3]45( r 4(*K(*K(*)4(* / ))tf> (2) 

where «[23]45 = ^(^2345 — ^3245), ^23 = ^23 + U 23 {t) and the notation U in the GF F 3 i/(tt'\U) means that it satisfies 
the equation of motion in the external field U. The boundary conditions for GFs are based on cyclic permutations of 
operators under sign of trace and, therefore, are not changed by the non-orthogonality. 
The Dyson equation for the full GF, 

{6 13 id t - or 2 1 / d*i[(&2 3 (*) 

-v%)6{t-t x )-V 2 z{tM\U)]}F 5 v{hJ\U) 

= tS(t-t')0- 1 ), (3) 
can be written if we write the double-electron GF via a functional derivative 

(T4(t)a4(t)a 5 (t)al,(t'))u = [F 34 (t,t+\U) 

d ^ F 6V (t,1?\U), (4) 



and define the self-energy as 



SU 3i (t+) 



Z 2V {t,t'\U) = { V[2 3]45 [F M (t,t + \U)+ 

F^txlU^F-titx^lU). (5) 



SU 34 (t+) 

Here, the inverse of the GF is defined by the equations 

F 56 (t,h\ U)F 6 ~} (h,t'\U) = F^ 1 (t,ti\U)F &v {t x ,1?\U) 

= 8(t-t')S lv . (6) 

The Hartree correction appears if we neglect the functional derivative in the self-energy 

Xg,(t,t 1 \U)=v [23]il ,F 34 (t,t+\U), (7) 

and the Fock correction comes from the derivative SF° /SU : 

E&,(t,ti|E0 = ~v [23]V5 F 53 (t,t+\U). (8) 

As we see, the first-order corrections to the self-energy do not contain the overlap matrix. The next corrections does 
not contain it cither. We will also need the RPA-screened Coulomb interaction. The difference between the standard 
procedure for Fermions and the one in this representation consists only in the necessity to write the representation 
indices. Introducing the effective field 



1 The external field U23 should not be confused with the Hubbard U 
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£#/(*) = EM*) + E&(t) (9) 
and S' = X — S H we can rewrite £' in terms of U e ^ as follows: 

y' f^'im v SU 7l f &) SF 5e (t, tl \U) F -i ( . ,n m 

<^78 (*2j 

= few , t + )J*9 (t, fa 1 CO 5F ltn![ W) ■ (io) 



The effective interaction, u, is 



5 2 785(t 2 ,t + ) = W [23]45-^7^7J = « [23]45 ^34 (*2, (H) 

Using the zero approximation for the vertex (F~ x — > T^T 1 ) 

5U 7 & J (t 2 ) 

we obtain the desired expression for self-energy in RPA: 

£^ A (M'|[/) - ^271'5(i\t + )^57(M'|t/)- (13) 

The inverse of the dielectric permeability, e -1 , entering the effective interaction is determined by the equation: 

e^, M {t',t + ) = 8{t+ -t')8 37 8 V i 

+V [73 ]4V F 35 (t,t 5 \ U)F 6i (t 5 , t+ | C/)£ 56 1 34 (t',t+). (14) 

Putting the external fields {7 = and making a Fourier transformation with respect to time we find the equation in 
terms of Matsubara frequencies (below we use ioj for Fermionic and iQ for Bosonic frequencies) 

ejy 34 {in) = S 37 S V4 + v [73]4 vU { °l 64 (in)e^l 34 (in), (15) 

where 

U^ M (iQ) = Tj2F 35 {iuj)F M (iLU + tfi). (16) 

UJ 

Then, the self-energy (U = 0) is 

S^ A (iw) = -v [23]45 T^2e 7 ^ 34 (iu - twi)F 6 7(iwi). (17) 

Let us write it as sum of "static" and "dynamic" parts: 

s^ A M = + s^ 21 ,m, (is) 

where 

^ = -«[23]46e s - t } 7 l' ,34 T E F 57(^l) = «^ , (19) 
S dyn A 21'M = -V [23] 4 5 TY,[ £ n',34( iw ~ iuJ l) 

Ull 

-^t;71',34]^57(^l). (20) 
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Here we separated the frequency-independent part of the dielectric permeability £^.72/ 34 . Its contribution to the 

self-energy, S^ p ^,, which is the frequency-independent part of the screened exchange interaction, is denoted v^y . 
The formulas given above have to be compared with the corresponding expressions in SCPT. We will see below that 
in SCPT /- and non-/-clcctrons are renormalized differently, for this reason we will need these formulas in terms of /- 
and non-/ indices (c) separately. Each of the corrections, say, E^,,S^ A contains four terms since each index takes 
two values: 1 = (c, /), where c = (jL) and / = (j/i) 

= V[jLJ 3 L 3 ]j i L i ,j'L'Fj 3 L 3 ,j 4 L 4 (tt + \U) (21) 
+ V yL,j 3ti3 ]j i L i ,rL'F 33 ^ 3 , u LA U+ \U) (22) 
+ v ijL, j3 L 3 ]3 4tM ,j>L>F j3 L 3 ,j itJ4 ( U+ \U) (23) 

JrV [jL ,j 3 fj. 3 ]j4fj.ij' L'Fj 3 ti 3 ,j4.tii (tt + \U), (24) 

= -{vjLJrLrJ'L', hL 5 (t't + )F hL5drL7 (tt'\U) (25) 
+v jL , j7 » 7 ,fL>, j5 L 5 (t't+)F <hL5th ^(U'\U) (26) 
+v jL , j7 L r ,f L >, h»At't + )F <hmdrLr (tt'\U) (27) 

+VjL,j r UT,j'L', j 5t i £ ,( t ' t+ )F j5 ^ 5 3 7 ^ 7 {tt'\U)}, (28) 

and the same for the other contributions £^ j, L ,{tt + \U), S^f A , L , (tt'|C/), etc. 

Thus, we can conclude, that within WCPT the overlap matrixes enter only the equations via the definition of the 
zero GF. As we shall soon see, the situation is more complex within SCPT, since in this case, the overlap matrixes 
enter in a form combining the inverse of the overlap matrix with many-electron population numbers. In the case when 
we take into account only the static part of the screening, the formula for the total energy can still be written in a 
simple form. Let us write it on the basis of Eq.(|36|). Let us denote the Hartree part of the Hamiltonian as 

ftf 2 = (l|fc ff (r)|2) = {l\^ + v ne {v)+v H {v)\2). (29) 

Then, representing the overlap matrix in the Cholesky form, O12 = Zi 7 Z 7 2, and introducing the eigenvectors, u", 
and eigenstates by the equation 

E Z nl + ■ < = (30) 

125 

We see that the solution to the equation for the GF in this approximation, 

{0 12 uj - [hg + v[f]}F 2V (w) = 6u> , (31) 

can be written in the following form 

= ((«il4»L" et) = E Z ~^ <U :^ . (32) 
u> — e v + 10 

Then, the expectation value (a| a i) is given by 

(4a 1 )=Y,Z^uy(sl<Z- 1 1 . (33) 

aryu 

The charge density p(r) can then be written as follows 

p(x) = -i lim V^WtVe^^Hl^^x') 

x'— >x * — ' z — * 

12 LJ 

12 



-* lim E^i(^E e iu ° F° FT (iu>)]4>l{x') 

x'^-x * — 4 * — 4 
12 uj 

J2^')M^W) DFT - (35) 
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Further, we can use the Migdal-Galitskii expression for the total energy, written in terms of, say, retarded Fermion 
Green functions, 

E tot = \J &uj[0 12 uj + (h° 12 - vg)]f(oj)(~) 

xImF 2 [ [ et \uo + i8). (36) 

In order to avoid possible confusion we note that F\ 2 (iio) are the coefficients of expansion of the GF in a series of the 
functions ^^(r), 

F(r,t;r',t') = -(T^(r,t)\^(r,t')) = 
i 

Y^M^F^W^x 1 ) (37) 

12 

and they do not coincide with the matrix elements F\ 2 {t,t') = J drdr'^^r)^^ 6 ^ (r, t; r', t')<p 2 (r) calculated on this 
functions. 

Inserting the imaginary part of the GF, Eq.(|32|), into the expression for total energy, Eq.(|3^), we find within the 
static RPA 

EUt = \ E Z^[Ox2e" + (1|^ + v ne (r)\2)} 

x f[?)v%Z-l (38) 

= \[{l\f- +v ne (r)+v H (r)\2)-v{f 
2 2m 

+ (l||^+^ e (r)|2)]<a t 1 a 2 ) (39) 
dr(^(r)[^+ Une (r)]^(r)) 



1 

+ 2 



(40) 



as expected. As follows from the comparison of the standard Kohn-Sham expression for total energy used in DFT 
calculations, the term v^ 2 \a\a 2 ) should be identified with the term J p(r)&-,(r) in DFT, i.e. with the contribution to 
the total energy from the screened exchange potential. According to Kotanio this expression in a static approximation 
(u> = 0) should give (at least, for normal metals) values close to what is usually obtained within standard LDA DFT 
calculations. 

One can write the expression for E to t also for the case of frequency^dependent screening. Then the self-energy also 
depends on frequency. For this case one can apply the idea of MigdalEj to introduce energy dependent eigenfunctions 
in order to diagonalize the GF 

E + ^£? 2 PA H]Z 2 7 • uJfM = (41) 

125 

Then, the Fermion GF is 

F (ret *> - tin l„t\\(ret) _ ^"SH"™ 'M^ai* .. . 

a~yv v ' 

The spectrum of the single-particle excitations is then determined by the poles of the GF 

e v = E v (e v ), (43) 

where 

E»{e») = u*"(s»)Zll[h? 2 + Re^\e»)]Z-lu^). (44) 



G 



In turn, the total energy acquires the form 



&tot = 2 Z 2 7 u 7 ( £ J 



I27Q 



xf(e")uZ(e")Z-l. (45) 

Thus, the eigenvectors ^(e") are needed only at uj — e u and the renormalization factor 1 — ( ^g^^ j takes 

into account the correlation induced deviation of the distribution function from a purely Fermionic one. 

At last, we have to note that in the case when one takes into account also the decay of the quasi-particles, the total 
energy can be calculated directly from Eq.(|36|). 

The facts which we are going to exploit further are: 1) the Galitskii-Migdal formula, Eq.(^6|), gives the total energy 
in the approximation chosen for Fermionic GFs; 2) the Fermionic GFs can be expressed in the form of a linear 
combination of the GFs for the Hubbard operators; 3) the lower intra atomic Fermi-like transitions give only core-like 
contribution to the energy. The latter statement will be considered in the next section. In the region of strong intra 
atomic coupling, the self-energies Sj^j'^', Ej^tj'z' for Fermionic GFs also can be found for many cases by making 
use of the SCPT where a non-perturbative calculation should be used for the interaction between electrons in each 
ion, but perturbative ones with respect to the interactions connecting different ions. How to calculate the GFs for 
Hubbard operators is discussed in next subsection and the details are given in Appendix |A|. 



B. SCPT and the Approximation of Renormalized Fermions 



If the Hubbard intra atomic repulsion is strong enough the perturbation theory from the atomic limit generates the 
following physical picture. The partially filled d- (or /-) shell is separated by a large energy gap to two subs hells. One 
of them can be viewed as core-like_states, while the other one can, in a certain approximation, be described in terms of 
effective, renormalized, FermionsE.3. These quasi-Fermions behave very similar to the Fermions in the weak-coupling 
regime. However, this simple picture arises only within the lowest approximation where the upper Hubbard sub bands 
are similar to the corresponding sub band of the so-called Hubbard-I approximation in the Hubbard models ^| and 
the decay of the quasi-particles is only due to normal scattering of these quasi Fermions caused by Coulomb intersite 
interaction. The scattering on collective excitations like spin waves, etc., as well as the scattering caused by .kinematic 
interactions, is not taken into account. In this case the Fermi-liquid type of decay (~ T 2 ) still holdsE^O at low 
temperature and excitation energy. Our aim will be to show that, in the ab initio calculations, for materials where 
the strong correlations are not developed very much due to the suppression by the a strong localization of a part of 
a shell, one can use the standard expression for the exchange-correlation potential even in the case of strong electron 
correlations (SEC) with the difference that a modified expression for the charge density should be used and, also, 
hopping, mixing and overlap matrixes should be slightly renormalized. In general case the renormalization constants, 
however, cannot be found without the system of equations for the GFs derived within SCPT. After some simplifying 
assumptions only one constant enters the equations and, therefore, only one equation should be added (see next 
paper). 

The idea of the derivation performed consists of the following. The density functional theory is based on an self- 
consistent solution of the Kohn-Sham equations for the charge density and potential. The key question is, of course, 
from where does one take the analytical form for this potential. Usually it is taken from the theary of the homogeneous 
electron gas and phenomenologically extended to the case of non-homogeneous systems. ShamQ, using the Hohenberg 
and Kohn theoremlij, has suggested a theory which connects the exact exchange-correlation potential with the exact 
self-energy. Although it is far from obvktus that the same holds for the case of an approximate self-energy, we will 
use this assumption. Sham and SchluterEj managed to solve this equation for the case of semiconductors within the 
approximation of spherical charge spheres and found that the analytical form in this approximation is the same as 
in the case of an phenomenological extension of the theory of homogeneous gas. At last, Kotanitj has shown that 
the results of the standard DFT-LDA calculations for the total energy are well reproduced within the static random 
phase approximation (RPA). These facts can be used as follows. If we manage to 



2 Although Hubbard did not consider in his works the inter-site Coulomb interactions, the approximation we use has the same 
underlying physics and, therefore, we will use the same name "Hubbard-I" for this approximation. 



7 



1. establish a one-to-one mapping between a subset of the Feynman graphs for the self-energy within SCPT and 
WCPT (at least, in RPA) and 

2. show that the Sham equation for the exchange-correlation potential within this subset of graphs has the same 
form as its analogue in WCPT and 

3. derive the expression for the charge density pp{r) within the approximation in the SCPT, which is the analogue 
of the RPA in WCPT, 

then, the analytical form for the exchange-correlation potential in this approximation should be the same as in the 
standard case, but the usual expression for the charge density p(r) should be changed to pp{r): v xc [p(r)] — > v xc [pp(r)], 
as well as the matrix elements of different interactions should be renormalized as it is required by the mapping. Since 
the LDA for the exchange-correlation potential works well for the case of delocalized electrons, we, thus, have to show 
that any graph in WCPT describing the Coulomb interaction between the electrons which belongs to one site and 
the other either to another site, or to the delocalized states, has its counter-part in the SCPT. As follows from the 
formulas for total energy and the Sham equation, both require knowledge of the full Fermionic GFs. Therefore, we 
have to derive the equations for the GFs, G, constructed in A-operators and, then, express the Fermionic GF, F, in 
terms of GFs G. 

In the SCPT the /-electron operator can be expressed in terms of Fermion-like intra-atomic excitation. In the 
approximation Hubbard-I (HIA) the latter can be described in terms of renormalized electrons. The mapping looks 
as follows 

y Coul . ^ = v Coul . ja^a ^ (^CW^/p^ . f*(X a /VP°~) 

= v c a ™ 1 ■ f a , (46) 

where f a is the operator of an effective Fermion, and P a is the spectral weight. Formally, this approximation can be 
obtained by a simple exchange of the commutation relations 

{A a , X h } = efz* = ef(Zt) + [efz* - ef{Z*)] 

~ef(Zt), (47) 
{ Cl ,xl}^O^f%ef(Zt). (48) 

Within the orbital representation where e^ 6 = S ab S(^ — [T, T]), i.e. it can only take the values 1,0, and the matrix 
elements f% = 1,0,-1. The expectation value is non-zero only when it is one of the population numbers, (Z^) = 
^,[r,r](^ r ) = 5^ >T] N r ; then 

ef(Zi) = ef(S arx] + S Ujtl] )S(a b)P a 

= 5(a - b)[N T + (49) 
for the transition a = [7, T] and, therefore, 

{ Cl , X\} ~ 0£ 2/ fiP?6{a - b). (50) 

Then, the A-operator can be written as X = (X a /V P a )^/~P^ — ] a \J P a . Then, the interaction is renormalized by 
the square root of P, which is brought by the X- (or, due to non-orthogonality of the basis set, c-) operator into the 
equation of motion for the GF and in the lowest approximations plays the role of the spectral weight of the pole, 
corresponding to the transition a. If the fluctuations described by the GFs, (T[efZ^(t) - ef(Z^(t))][efZ^(t') - 

e^ b (Z^(t'))]), are neglected, the perturbation theory contains only Fermi-like Feynman graphs. The latter means 
that any Feynman graph obtained within a weak-coupling expansion has its counter-partner in the strong-coupling 
perturbation theory. Attention has to be paid to the fact that the number of renormalized Fermions does not coincide 
with the number of ordinary Fermions: the number of the latter is equal to the number of orbitals involved, while 
the number of the former is equal to the number of Fermi-like transitions involved. Although it is not difficult in 
principle to answer the question, how to derive this approximation and what remains beyond this approximation, it 
requires quite lengthy calculations. 

Different ways may be used for the calculation of the Fermionic GFs within the SCPT: 1) directly, cither via a chain 

of equations of motion, until the chain becomes closed on one site via the operator relations / 2 = 0, 

equivalently, in terms of chain of high-order GFs; 2) using the representation of Hubbard operators, which diagonalizc 

the on-site interactions, and the diagram technique for them. Then we have to find the Fermion-like GFs via the 
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Hubbard A-operators. Here we use the second way, since the SCPT for the GFs are needed and some approximations 
are considered in a paper by Sandalov et aZ.EJ Let us inspect how the first-order corrections to GFs generated by the 
same types of the matrix elements of the Coulomb interaction look like (one should remember that the inverse of the 
overlap matrix is not included into the definition of self-energy). 

The term of the Hamiltonian describing the Coulomb interaction between electrons can be written in the jL- 
representation in the form 

Hcom = ^2345 [4 + (Ar 2 XrH + (flf 3 X^} 

x [ci + frxHln + fFxs*]. (51) 

The interactions with the core electrons are included into the definition of the A-operators here. Now we have to 
consider all these sixteen terms in order to take into account single-site correlations via multiplication rules for the 
Hubbard operators. It is convenient to consider the terms which contain different number of /-operators separately. 
Thus, we form classes according to the number of /-operators contained in the interaction term. In order to reduce the 
number of terms in the equations of motion we collect similar terms by making use of the commutation relations and 
a proper symmetrization of the interaction. We also use the short-hand notations for indices: for conduction electrons 
l c = (ji,Li) = (ji,h ^ 3,mu,si = 1/2, oi) and for /-electrons 1/ = (ji,/ii) = (h,h = 3, mu,si = 1/2, 01) , where 
it does not lead to a confusion. Also, it is convenient to introduce the group constants which are the coefficients in 
the multiplication rules for A-operators: 

X" • X) = nfzl, X? • Z] = nfX*;, etc. (52) 

Obviously, 



ef = nf + nf (53) 

and so on. 

The /-operator in ^-representation is given by a sum of A-operators, and we are not able to find the Fermionic 
GF FVt) = {Tfi{ri)fl,(T')) without finding the GFs G Xc , G cX , G xx , involving Hubbard A-operators. Therefore, 
we need the equations for these GFs. For the Hamiltonian which includes hopping and mixing interactions (Hubbard- 
Anderson model) these equations are given in a paper by Sandalov et aZ.Ell Therefore, we have to add the terms 
generated by Coulomb interaction. If we write the equations in terms of functional derivatives with respect to 
the external field, all mixed terms arise when we iterate these equations. However, as seen from the form of the 
contributions from the Coulomb interactions to the equations of motion for the operators c and X (see Appendix ^), 
the fields which are introduced in Ref.Eil, are not sufficient. Now instead of the S^xt-matrix, describing external fields 
we have to use the following S cx t 




S cxt (-i0,O) = exp{ -i I -//£,,,(/) } . (-vl.) 
where 

jL.j'a' 

+ J2 xf(t)u x : fL ,(t) CfL ,(t) 

ja,j'L' 

+ Xf(t)U x a x a ,(t)Xp(t)(l-S jf ). (55) 

jaj'L' 

As seen, all terms except the one proportional to TJ% renormalize, the c-c-hopping of the conduction electrons, mixing 
and /-/-hopping: 



H cc H cc = H cc + JJCC 

W = W + U cX ', 

t^t = t + U cc , (56) 
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therefore, we consider that all these external fields are included into hopping and mixing matrix elements below. 
These fields generate the terms already considered in a paper by Sandalov et al.c3, in perturbation theory, and the 
problem has a non-linear nature from the very beginning. Now, using the the contributions into the equations of 
motion from the Coulomb interaction, given in Appendix [X], we will express the new complex GFs in terms of the 
simple ones. It is convenient to consider these contributions class by class; the number of the class corresponds to the 
number of /-operators in the term of Coulomb interaction. Before doing this it useful to note that contrary to the 
weak-coupling expansion the representation which we use here is not unique. However, after one is chosen, it generates 
the system of graphs of the diagram technique which corresponds to this choice of the closed form of the equations 
for the GFs since we obtain the graphs by means of iterations of these equations. As a consequence, special care is 
required in comparison with some other techniques: one should not expect graph-to-graph correspondence between 
different diagram expansions for the GFs. Our convention here is the following: in the first step we separate the first 
Bose-like operator which stands in the left-time product and then, if necessary, make it again. This sequence simplifies 
the calculation of the necessary time limits. The analogue of the Hubbard-I approximation (HIA) is obtained if one 
fully neglects all the contributions coming from the functional derivatives. 
The approximation of " renormalized Fermions" can be introduced if 

a) we will show that, from a physical point of view, the system of correlated electrons has the Fermion-like exci- 
tations, which may differ very little from the ones in an weakly interacting Fermions gas, at least in some region of 
parameters (strength of interactions, dimensionality, temperature, etc.); 

b) from the point of view of mathematics we will show that even in the strong-coupling regime a subsequence of 
Feynman graphs in SCPT exists which describes these excitations. Then, the remaining graphs will determine the 
region of parameters where this approximation is valid. 

We will now search for a solution to the equations of motion for the GFs in the form 

G = VP, (57) 

where the bold letters denote matrixes with respect to all indices and times as well as their product implies matrix 
multiplication with respect to all indices. Here the GF is 

G aP = \{Tr, a {t)^H)) u% rj a (t) = c jL (t),X?(t). (58) 

The cc-component here G\\ — G^ cc ' = is a Fermion GF. In terms of the operators 77, the interaction can be 

written as 

1 - ,t„t. 



Vcoul = 2 W 2345% ? /3 7 74 7 75> ( 59 ) 



with 



v 2345 - V j 2 L 2 ,j 3 L3,j4L i ,j 5 L B > 
—c fee 



(ft \a 3 . 



2345 jzLatjaHaifaLtthLs "^ 3 y 

S8&° = -nL^u^LMLriU.r, etc. (60) 



The external fields in Eq.(p54|) and Eq.(|55|) can also be written shortly as 773D34774. Of course, the operators belonging to 



the same site should be multiplied according to the rules for the A-operators, Eq. (p2[) . We shall consider these terms, 
later. Since the anticommutators {c, X' }, {A, X'} give an operator, not c-number, we have to write all equations for 
a general case {771,772} = e\ 2 Z 3b , where Z is Bose-like operator and, therefore, Vcoul gives the following contribution 
into the equations of motion 

[Vi, %oui] = w [23]45 £3^ 3b ?7^77 4 77 5 . (61) 

In the case when 771 = c\ and r\\ = c\, the operator Z 3b — 1 and £3^ = \ , In the equation for the GF the term 
generated by [771, Vcoul] can be written as follows 
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y [23]45 £ 3i, 



(Z 3 *(t+) +1- 



SU 3 »(t+ 

x{Tvt(t)V4(t)v 5 (t)r)Ut'))u 
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u [23]45£ 36 



(Z^{t++))+l 



(T V l(t + )r li (t + ))+ l ] 



SU s »{t++)_ 

GH{t,t'\U). (62) 



8U 34 (t+) 
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Here t ++ and t + denote the limits limt 1 _,t+o li m t 2 ^ti+o respectively. Recall that [FQ C (t, t' 1 1/) ~ 1 ] 51' does not depend 
on the fields U 3b and, therefore, SF^,(t,t'\U)/8U 3b (t") can give non-zero contributions only via the mixed GFs 
(TcX*). Therefore, we can write the equation for the GF in the following form 



x [V23 + U2 3 }G 3V 



[Z ib )+i 



5U 4 



-V „f\ 2 



su. 



4*, 



J [23]45 



-G 43 (t+ ,t 3 ) + i 



5U 34 (t+) 



"51' J 



(63) 



where for brevity we write the time argument only in those places where a confusion may arise. The mixing and 
hopping are denoted by V following the definitions in I. We see that there is an essential difference with the standard 
WCPT: here we have to calculate the functional derivatives twice, while in WCPT only second square bracket is 
present. Now let us, in order to give the reader a feeling of the SCPT, make the first iterations and give both the 
analytical and graphical expressions for the corrections to the GF. In zero order of SCPT G^, G"j, = T>\ V iPy,y. 
Let us now insert G\\i =>■ G°j, into the right-hand side. We have the terms where S/SU = 0, 



5G { V=Vl v 4l\z^)[V- 23 + U n ]Gl v 

-Vl v ,eXlHz ib )v m ^G A3 {ttM)Gl v . 



(64) 



Since e 1 J' 2 (Z ib ) 



P 1 2 and G = VP, we have simply 



5G^ = GUV2 3 + U- 23 ]G% 



1"2 



v {2 3 ]i5G%{tt ,t 3 )G 5l , 



(65) 



We denote the pseudolocator D by a solid line with arrow in the Feynman graphs, the end-factor P with an open circle, 
for the mixing- hopping plus external field, V + U, a wavy line is used and for the Coulomb interaction we use a dashed 
line. Below, we will also need a notation for the Bose-like correlation function K lb2b (t,t') — 5(Z lb (t)}/ SU-jAt'); we 
will use the curly line for it. The graphs Fig. [TJ and Fig. || correspond to the analytical expressions in Eq.(|65|). Since 
SG^/ SU 3i = and this derivative removes one of the interactions [V23 + ^23] m anv expression, in order to obtain the 

other graphs of first order with respect to Coulomb interaction, we have to insert G§y ^> Sy^G^y = G^^V^ + U^G^, 

into the generating equation. Then, the term with the derivative 6[5y Ggj/]/ 5U 3i produces the exchange graph (see 
Fig. |): 



Gi2 V [23\i5G%G% 

Next, we find four terms generated by the derivative 5/8U± b : 



-1"2 = 



'[23]45' C> 46 (*3 ' ^6^67 \4 b G% (*7 , h)G^ j, 



ffl 1 "2- -r-,0 63 K r 

u ll" £ 4: b u [23]45^46 e 6 i , A 6 i) 4 i) (j 51' 

V { l v ,e\ h 2 v [n]A5 G%{t^ ,t 3 )Vl & Tl 7Ab G a n , 
^v4i% 3]A5 G%(t+M)Vl & K 6b , b e 



61' 



(66) 

(67) 
(68) 
(69) 
(70) 



The corresponding graphs are shown in the Figures ||, |^, ^ and R respectively. Continuing these iterations we find 
that there are a subsequence given by the graphs in the Figures g, [lO] and 11 , which exactly have the structure 
of the WCPT: the only difference is that instead of Fermion GFs, F, in WCPT here we are dealing with the GFs, 
G. Obviously, each one of these GFs, G, can be dressed with hopping and mixing, as shown in Fig. |l3[ if the 
corresponding transition is in the energy region where these interactions are not equal to zero (see next sections). 

Then, these GFs describe delocalized transitions, or, in other words, Hubbard type of bands. The graphs [|, |[ ^ 
and ^ do not appear in WCPT for the single-electron GFs and describe contributions from kinematic interactions. 
Let us now return to the definition , Eq. (pOj) , of the Coulomb matrix elements, V[2 3 ]4Fig.5j an d consider, for example, 
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the graph ||. In the case when one of the indices, in Eq.(|67|), say, 5, describes a c-electron, there are no additional 
factors, f as . If the index 5 correspond to an /-electron, then there is a factor f as . Such factors are automatically 
provided by the matrix elements ^[23)45, Eq.(|60|), for all inner lines, but not for the external ends. Therefore, if we 
multiply these type of graphs by the external factors, say, for c, and use the expansion of the Fermi operator in terms 
of the Hubbard operators, we find that 

f ai^fO - r<o_r<o_ /ft \ a i> — pO p°_ p°_ /7i\ 

J vi ^12 "[23]45 53 ^41' Uiv J ~ r 12 y [23]45" r 53^41' ' \' L > 

i.e. these graphs describe the standard Coulomb scattering of electrons, but the Fermion GFs should be found within 
SCPT. Besides, as seen from these iterations, any application of the derivative of the sort 5/8U 3b (t") to the expression 
containing at least one GF G(t,t'\U) leads to the graph describing the effect of kinematic interaction and, therefore, 
destroys the one-to-one correspondence between the WCPT and SCPT in the ARF series. Thus, in order to obtain 
the equation for the GF in ARF we have to neglect in the first bracket of the Eq.(||) SG(t, t'\U)/5U 3 »(t"). Then, the 
interaction is renormalized by the expectation value of the Bose-like operator 

^2345 - V 23 4 5 4 2 b (Z 3b (t ++ )) = V 2345 P 12 . (72) 

This condition removes not only the kinematic interactions, but also any of the graphs containing the correlation 
functions K (see the figures | and @). It is worth to note that some of these correlators describe spin waves in the 
magnetic ordered media, therefore, this type of scattering of carriers is beyond the ARF. In spirit the ARF corresponds 
to the Hubbard-I approximation (in our definition). 

This remarkable similarity can be achieved in the level of equations of motion for the GFs. Although this equation 
is approximate it is still a functional equation, therefore, the similarity holds in any orders of the PT, which can be 
generated by iterations of these equations. Therefore, the analogue of the RPA in WCPT can be constructed too. 
This is obvious from the graphs in the figures |ll| and [l^ and further graphs containing an increasing number of loops. 
In matrix notation the definition of this approximation arises from the "reduced" differentiation of the GF. Namely, 
the higher correlation functions arising in the equation of motion are expressed in terms of functional derivatives of the 
GFs with respect to the bosonic fields which are needed to reproduce these higher order GFs. Thus, we will call the 
described approximation the "approximation of renormalized Fermions" . The factor P, is specific for each transition 
constant, renormalizing the hopping, Coulomb and mixing interactions. One can expect from this comparison that 
those abnormal features of some compounds which are different from the normal Fermi liquid have to be described 
by some of the remaining graphs which are not included into the set defined above as ARF. The analogy can be 
continued, for example, one can introduce the self energy in SCPT within the ARF in the form 

E SCPT = [vP((^r,)+i-^)G]G-\ (73) 
which is in complete analogue to the standard definition used in the WCPT, 

^CPT = [v{{p)+i _S_ )F]F -l_ (?4) 

Thus, we have established a one-to-one correspondence between the perturbation theory for smg/e-electron GFs, F, 
(weak-coupling expansion) and a sub-series in the perturbation theory for the ?7ianj/-electron GFs, G. This is the 
central result of the present paper. One of our targets formulated in the introduction was to build a bridge between 
the standard DFT-LDA calculations and SCPT and to find an extension of LDA which is reasonably simple and 
takes into account some of the features of SCE. The form of the exchange-correlation potential w^ 6 ' 9 ' which is used in 
DFT-LDA is derived from the theory of the homogeneous electron gas, i.e. far the case of fully delocalized electrons. 
Therefore, the potential v DFT , should be the solution to the Sham equationQ, 

= ^e i » 0+ {Fg FT (iu ) )[£ 23 (iu,) - vST]F 31 ,(iu;)}. (75) 

w 

is expected to have an analytical form close to v^ 9 ' only for the delocalized Bloch electrons. The connection between 
the matrix elements £25(10;) and the coefficients of expansion £34(1^) is given by the relation 



£25^) = / dxdx'(j> 2 {xi)Y;(x, x' , iuj)(f>l(x2) 

= 2 3£34M045, (76) 
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In turn, this means that the one-to-one correspondence, equations [74] and [73], found between the WCPT and SCPT 
series of graphs can be constructively used for the delocalized transitions. The recipe follows from comparison of 
these two equations: the following replacements should be made in the expression for the self-energy: F — > {f]ry) and 
the matrix element of Coulomb interaction v — > vP. Besides, the charge density should be rewritten in the form, 
where each /-operator is represented in terms of A-operators: / M = f^X a . As follows from I, and will be shown in 
details in III, the matrix elements of mixing interaction and overlap matrixes are also renormalized by the spectral 
weights. Let us now compare the equations for the dielectric permeability within WCPT, Eq.(p^), and SCPT. On 
one hand, according to our comparison of graphs, Eq.(|64])-(|7l|), and the more general observation, Eq.(|7^)-(|7^), in 
order to transform a graph of WCPT into the one of SCPT we have to replace each bare /-GF as follows 

0-2 

On the other hand, the low-energy transitions does almost not form bands at all. Let us write down the elementary 
loop (the graph beyond exchange) and separate the lower transitions in it . Each index in the loop 

n i°2 ) ,34(^) = Tj2 F i2(i^)F 34 (iuj + ifi) (78) 

f 

can take the values 1 = l c , ly, ly , i.e. the index of a conduction electron, j\L\, or one of values for an /-electron, 

j\v\ or j\V' x , where e denotes an empty (delocalized) and / stands for a filled (localized) orbitals. If we neglect, for 
a first step, the contributions from the Coulomb matrix elements with an odd number of f-operators^|, then the lower 
pole contributes only in the loops, which contain one GF of the //-type, i 7 i / 2 / or i 7 3 / 4 / - All other loops are standard 

ones with the only difference that the standard /-GFs, F\ 2 , when the f-electrons are treated as valent, however they 

acquire the factor \(f v J) a *\ 2 P a2 where the sum is taken over the upper transitions 02, and Eg should be replaced by 
A &2 . Since in the loop 

n ( % (til) ~ Y B, B 3 4 (k) f{A ^ } - f (e ** ~ ^ (79) 
/(A 5l ) is always equal to 1 and the Fermi function f(s% x — y) = only if e a kx > fi , the loops 11:°, f (iO) have 

1 jr 2 j. , 3 C 4 C 

large denominators in the region where the numerator is not equal to zero and, therefore, they are small compared to 
11^2 3c4c (if2). Here, the coefficients B transform the GFs from the orbital representation to the bands kX. Thus, as 
it is expected, the main contribution to the screening comes from the conduction electrons including the delocalized 
upper transitions. Within this approximation the formulas: Eq. (|l5|) for the dielectric permeability e -1 , Eq.(p"7j) for 
the self-energy vjRPA remams valid, however, with the replacement of the full GF by only the upper-pole part. There 
are, of course, also contributions from the electrostatic interactions between electrons of different sites, which are 
treated in a standard way with the difference that each transition enters with the its spectral weight. 

Let us turn now to the Galitskii-Migdal expression Eq.(|3^) for the total energy. In the last term of Eq.(|36|) we 
have a summation over the transitions. For the lower transitions Oj t j 2 — Sj 1 j 2 , tj 1 j 2 — 5j 1 j 2 e{ j and the local terms 
give a contribution to the energy from the core-like levels. The remaining expression, within the ARF, coincides with 
the contribution to the energy from delocalized electrons within the WCPT if we make the replacement in the /-GFs 
discussed above. Exactly the same arguments can be applied to the Sham equation, Eq. (|75|) , connecting the many- 
body GF and the self-energy with the DFT GF and the exchange-correlation potential. Therefore, we can conclude 
that the analytical form of the exchange-correlation potential used for the description of delocalized electrons can be 
appli ed also to the case of correlated electrons if we will use the expression for the charge density Eq. (|35| in the form 
( |100| ) (the derivation is given in next section). Note that in the limit when the lower spectral weight for the localized 
electrons P = 1 and A a2 — > 00 we are back to the standard model for the lanthanides (see the next section), since 
in this case we have n /-electrons localized and the upper transitions do not contribute to the energy. Then, the 
possibility to describe photo-electron spectroscopy disappear. When P < 1, we have to solve a system of equations for 
the spectral weights additionally to a standard DFT calculation and bring these two calculations to self-consistency. 
The remaining important question, the magnitudes of Ai and A2 will be discussed in last section. 



3 We can do it since even differential overlap between localized /-orbitals and delocalized c-ones is negligible 
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Since only a part of the transitions is delocalizing, some difference arises between electrons and quasi Fermionic 
excitations, which we call renormalized Fermions. First, the numbers of Fermions and renormalized Fermions are 
different, since the latter is determined by the number of Fermi-like intra-atomic transitions, involved in the physics 
under consideration. Second, the WCPT generates the graphs describing all intra-core and /-/interactions, relating 
to the same ion, while they are already taken into account via the algebraic construction of the many-electron state 
in the GFs, G. The intra-atomic Coulomb interactions are included into the WCPT theory separately; as will be seen 
later, in SCPT the transitions, describing deep intra-ionic transitions are naturally separated by a big energy gap 
from the delocalized electrons and actually the discussed mapping involves only the delocalized states, providing, (if 
necessary) a description of the many-electron structure of the ion including the Hund rules. Third In order to make 
the two series, WCPT and SCPT, fully analogous in the Coulomb interaction, Eq. (|so|) , and to derive the ARF, we 
had to treat all terms in the way as if they belonged to different sites. Therefore, we had to add the graphs, which 
treat the correction to the terms in Eq.(|60|) with a few /-operators belonging to the same site, and subtract those, 
which within the ARF, are treated incorrectly. This difference should be included in the error of the ARF. Therefore, 
the accuracy of this approximation (or an applicability to the physics in question) is determined not only by those 
graphs which describe the different manifestations of the kinematic interactions, but also by the difference between 
the correlated and decoupled single-ion terms in the ARF. In Appendix [A|, we have separated all Coulomb terms into 
classes containing different number of /(d)-operators (or core-operators): the terms belonging to the class do not 
contain them at all, the ones belonging to the class 1 contain one /-operators etc. Due to the large number of terms 
this consideration is quite lengthy. The terms generated by Coulomb interaction in the equations of motion for Green 
functions operators are given in Appendix |A|. 

In order to estimate the error, due to this difference between the correlated and decoupled single ion terms, we 
have to take all the graphs which contain some of the sites, or all of them, coinciding and subtract the graphs with 
the same value of the matrix element of Coulomb interaction but with the analytical expression for the set of the 
GFs corresponding to certain choice of links between sites. The sites are defined as linked if they belong to the 
same site. These different choices are considered in Appendix [X| The interaction of the density-density type is 
expected to give the main contribution, therefore, let us consider the graph 6 as a example. It takes into account 
the Hartree contributions from the terms of the Coulomb Hamiltonian: — W[ J - 2a2 ^ 3a3 ] J4a4 j 505 Aj 2 2 A° 3 3 A° 4 4 A°. 5 , 

v 2 = vi J2 a 2 .j 3 L 3 ]j 4 L 4 j 5 a & Xl 2 c] 3L3 c HLi Xl 5 and V 3 C = Vfo^tL^Ltj^c^c^Cj^CjiLi, where all the sites j { 
are different. There are also terms of the Vf type where the sites js = ji or/and j 2 = J5. In term s of Appe ndix [a| the 



contributions from these interactions belong to the class 4 and are described by equations (A24) to (A27). We have 



to put S/SU Z = in these terms. The term v jaSadaiiajaaAdBa6 X%X%X%X% = V^Masj^mK^Xg^Xfi 
generates the following contributions to the self-energy T,j a j 5a5 

r,. - . - . \x.. K - S 3 a 4 f aa 2 1 /yi \ 

-5 jj3 K? a *e?{Xf*X>}], (80) 

whereas the term which has been added for inclusion into the ARF and, therefore, subtracted, is 

U J2 02,^303^304,^505 

-5 jh e£*{zf){X%X?)]. (81) 

The difference is 

Vj2 0a J3a3J3a4J505%3 [— ^ 3Q4 £^ + e^ a3 S a4b (Zp] (Xj^Xj). (82) 

It is convenient to use the equality 

K f a *ef = £ f 3 4 a4 , (83) 

which follows from the commutation relations between Hubbard operators. Then, the difference between the terms is 
proportional to 

{-(Zl)S ail ' + /s|"*)ef s (84) 

for diagonal £ = [T, T] this is cx (1 — P ai ). Now we will use an information about the spectral weights which is 
considered in details in the next sections. Both, localized and delocalized /-orbitals have upper (A" ~ Ai) and lower 
(A a ~ A 2 ) transitions . For localized orbitals and for lower transitions a 2 or/and b in Eq.(B2), the expectation values, 
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(Xj 2 Xj), are small (<C 1) due to the suppressed hopping and mixing in the region of energies at energies u> ~ A ai 
and big energy gap ep — A Q1 ; the latter suppresses admixture across the gap. For the upper transitions the localized 
orbitals themselves have in this region of energies small spectral weight and again (Xj 2 Xn <C 1. The weight P ai is 
close to unity for the delocalized orbitals near the upper transitions (A2), therefore, the small constant decreasing 
the difference in Eq.(|82|) is oc (1 — P ai ). 

It is also worth to discuss a possible confusion: it may seem that these speculations are wrong just for the simple 
reason that the equation for the GF in the WCPT is exact while both representations (single-electron Fermionic 
and many-electron) are equivalent; therefore, WCPT must generate some graphs which should not appear in the 
SCPT for an approximate equation for the GF in SCPT. Indeed, all this consideration might be performed within 
the WCPT, however, in this case one has to consider the theory, not for the single-electron GFs, F, in WCPT but, 
for the many-electron ones, Gp(x\, x%, x$, x n ), (here Xi — (r,-,tj) ) which correspond to the functions Gx- The 
functions Gp(x\, X2, X3, x n ; Rf) describe the relative motion of all n electrons in an ion and, therefore, contain 
much more information than we actually need. This means that to find them require much more effort than for the 
simplified functions Gx, which are constructed on many-electron operators, X(t), belonging to the same time, t. Of 
course, the same argument shows the insufficiency of the approach discussed in the cases when a knowledge of the 
dynamics of each electron is essential. This situation arises when the magnitude of the perturbation is so large, that 
the complexes are destroyed completely. In Appendix [b] we consider an simple example, illustrating the difference in 
the descriptions already in the level of the ARF. Now we shall consider how the standard model for the lanthanides, 
widely used in ab initio methods, arises in our approach. 

III. STANDARD MODEL FOR THE LANTHANIDES AND ITS EXTENSION 

It is instructive to consider the limiting case corresponding to the so-called standard model for the rare earths first, 
where the /-electrons are assumed to be completely localized and atomic like. We will show that this model sometimes 
needs to be modified, in order to explain all experimental data of the lanthanides. Let us first have a look at the 
structure of the zero Fermion GF in the atomic limit, which is given by 

F(«*)M = (T/ v /t) iw = ^l^ a 



E 



iu - A a 

|(r„|/„|r„ +1 )| 2 (7Vr„ + iVr„ +1 ; 



iu- {Ev n+X -E Vn ) 



(85) 



The differences of the energies of the ion states Ej- n+1 — E-p n ~ Er n — Ej* n _ 1 ~ U, i.e. of the order of value of Hubbard 
U; here n is the number of localized electrons of the element in question in a given compound. Due to the large value 
of U the transition Er n — Er n _ 1 = Ar„_ 1 r„ in many cases (particularly, in rare earths) is much below the bottom of 
the valence bands, while the transition Er n+1 — £r„ = Ar n+1 r„ is above the Fermi level (as seen from spectroscopical 
data, in many rare earths it is not very high above ep). The strength of the hopping and mixing depend on the energy 
of the states involved since the width of the inter-atomic barrier for electron tunneling is increasing with a decrease of 
the energy and, also, the intra-atomic states of an isolated ion belonging to different energies are orthogonal to each 
other; for the deep states of a non-isolated atoms this picture is changed only by a small amount^. On one hand, as 
discussed earlier (see Fig. 5, 6 in RefB), the mixing interaction and hopping gradually, but quite fast, are switching 
off when, say, the /- or d-level moves deeper to the energy region below the bottom of a band. In our language this 
corresponds to the lower transitions, Ar n _ 1 r„- Therefore, only the upper transitions, which are in the proximity of 
the Fermi level, can participate in the formation of bands. However, if the upper transition is much higher than ep, 
the bands formed are not occupied and do not contribute to the cohesive energy. From the physical point of view 
this picture is fully analogous to the standard model: in both cases the localized, occupied, /-bands contribute to the 
total energy as core states while the contribution to the cohesive energy in both cases is absent. 

The arguments given in the simple example in Appendix |b| can be easily extended to the general case of n orbitals. 
Let us show how this separation (or, in another language, "orbital polarization") arises. The immediate objection often 
arising is that from the point of view of simple models of strongly correlated electrons, say, Hubbard-like models, the 



4 It should be noted that this statement is representation dependent and relates to the methods where the intra-atomic wave 
functions in the basis set are taken as atomic-like. 
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switching off of mixing and hopping seems strange. Indeed, the transformation of the Hamiltonian from the Fermion 
representation to the representation of Hubbard operators is exact, 



HHubb = E ( £ °M A,y + tiDfLfj» + E u *i&» ( 86 ) 
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therefore, if one writes the same magnitude of hopping for all orbitals, = ty, its magnitude for the lower transitions 
coincides with the one for the upper transitions. Here we have written a truncated single-site Coulomb interaction in 
order to simplify the discussion; the Greek indices denote the orbitals: A = (I, mi, a), the vacuum state |0) includes 
all filled shells of ion, |A) = // A |0), |A/i) = /^/^JO), etc., and the orbital A belongs to an open shell. However, this 
is an oversimplification which makes the physics different from what happens in reality the hopping is always orbital 
dependent. Let us denote the occupied orbitals of the shell by Uf, and unoccupied ones by v e . In lowest order one 
expects that the Hubbard sub bands will be generated by the poles of the Fermionic GF. Indeed, the "band structure" 
in the Hubbard-I approximation for the model t^ = S Xu t^ is given by the equation: 
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- t v k = 0. (88) 



However, as discussed above, if the orbital ff is deep in the potential well, the hopping t^f becomes so small (expo- 
nentially) that the bandwidth becomes of the order of the width of the transition, T a (the width is due to interaction 
with phonons, electro-magnetic fluctuations etc.). In this region of parameters the coherence is destroyed and bands 
are not formed. Besides, the bare hopping entering this equation, actually, does not determine the real band struc- 
ture: the equation for it contains the self-energy which also should be found self-consistcntly (in both, WCPT and 
SCPT); the hopping itself in the Hubbard representation is treated as an interaction and is dressed either with the 
vertex corrections within field-theory methods or, via changes of the wave functions within ah initio band structure 
methods. Second, as follows from the observation that any element shows a valence which is either integer, or, for the 
compounds with intermediate valence, between to nearest integers in a given surrounding and pressure, only three 
groups of the transitions 
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(89) 

can contribute to the Fermionic GF. Here n is the number of localized electrons in the ion. In order to simplify the 
discussion we take |r„) in the simple form: |T„) = /!•••/« 1 0) , put A 5i = Aj, i = 1,2,3 and omit the widths of the 

levels ri ( . Then the matrix elements |/°f = K0|/„/„_ 1 .../ i ,_ 1 / I/+ i.../„_ 1 |M/ 1 t / 2 t -/ i ,-i/^+i-/^|0)| 2 = 1 and arc 
equal to zero if the orbital v is not present in the state |r„) and the equations for the spectrum become 
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Let us consider the combinations of the population numbers P a = P a ( r T ), which are the spectral weights for the non- 
perturbed Fermionic GFs. For simplicity we also take the sum rule for the population numbers in the form Nr = 1 
(recall that this sum rule is representation dependent; it follows from the commutation relation {fju,fj V } = 1, but 
this is valid only for the orthogonal basis set of wave functions; in the case of non-orthogonal set {fju,.f]\} = 
and the sum rule is also different). According to our assumption the ion has n electrons localized, i.e., the states 
^C^„, v n ft.ft 2 -~fl n \fy — |r^) are filled in the non-interacting ion. If the interactions which remove the degeneracy 
of these states are absent (not taken into account, or the temperature is higher than the corresponding splitting), 
then, there is no physical reason for the population numbers of these states to be different. Let us calculate the 
population numbers of the states involved for one chosen orbital and then use the symmetry between the orbitals. 
Let us choose the orbital v = 1. Then any of the transitions fi, v-i...v n -\) ~~ * l^i, ^■■■f-'n-i) for vi from the set 
v = 2, 3, N gives a non-zero matrix element I/" 1 1 2 ; therefore 

N k> = (n) = (93) 
and all other population numbers are zero. Then 

PSo = E \f% I' E £ [rx] N r = N$ n _ t + N° n = N Fn (94) 
ai r 

and 

Kl = E \fH I 2 E = ^ + N Fn+1 = ^r; . (95) 

oi r 

Therefore, the zero Fermion GF for the particular orbital v\ is 
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(96) 



The total spectral weight is equal to one as it should be. However, at u> = Ai we have the weight n/N < 1, while the 
weight of the low-energy transition from all iV orbitals is n, as it is expected within a single-electron picture. Thus, 
within the non-polarized picture one can expect from Eq.(|96j) that the bandwidth of the lower Hubbard sub bands 
will be proportional to n/N and the upper one to (N — n)/N. ^ This paramagnetic picture is, however, too rough. 
First, this speculation does not reflect the fact, established in the ab initio calculations, that if the energy A ai is below 
the bottom of the lowest conduction band, any matrix element of hopping or mixing, calculated within self-consistent 
scheme becomes negligible^. Second, the degeneracy is certainly overestimated: the first Hund rule is caused by the 
intra-atomic exchange integrals, but the splitting between different orbitals may be a few such integrals (say, the 
transitions between three-electron states and two electron ones cover the range from 2U — 2 J to 2U — 9 J, where J 
is exchange integral); the next-order Coulomb interactions, which are responsible for the second Hund rule, are also 
much greater than room temperature; at last, the spin-orbit coupling, for example, in rare earths provides a splitting 
of eV-order, between the multiplcts. Thus, if the atomic-like Hund rules are not fully destroyed we find in a zero 
approximation that only one of the population numbers ^Vp 2 ~ 1, and, therefore, due to the sum rule, N® = 1, 
all remaining = 0. The GFs 

pax 1 

uj- A ai uj - A ai 

Gi° 2 ) M = -^-«^^, (97) 
uj - A a2 uj - A a2 



5 Within the diagrammatic field-theory methods this problem in the case of degenerate states is not solved yet and we postpone 
this for a future work. 
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and the Fermion GF for the filled orbitals is 



(98) 



i.e. we are dealing with the situation, considered in the example in Appendix [b|. Thus, let us switch on the hopping 
and mixing. We have to distinguish the Hubbard sub bands originating from empty and filled orbitals since hopping 
for them are different. 

In the case of empty orbitals the hopping t„ e ^ since they have energies above Fermi energy, where the potential 
barrier between atoms either is absent or is small. However, as we have seen in the previous section, the dispersion of 
the lower Hubbard sub band is negligible due to the small weight of these orbitals in this region of energies. In the case 
of filled orbitals the dispersion in the lower Hubbard bands is also negligible since t Vj is small due to the large width 
and height of the energetic barriers between atoms; near the upper transition their spectral weight is small. Thus, 
the only dispersion which is physically important is the dispersion of the empty orbitals near the upper transitions. 
Above we have shown that the switching of the mixing and hopping leads to the appearance of non-zero spectral 
weights, above the Fermi level for the filled orbitals and for the empty orbitals in the region of lower transitions. But 
the dispersion in the latter region is negligible for both type of orbitals. Therefore, the role of delocalized electrons is 
played by delocalized upper transitions of the empty orbitals and we can use this fact in order to derive some recipe 
for how to correct the exchange-correlation potential. On the other hand, we can in the general formulas use the 
approximation t v ° ~ as well as the overlap matrix with neighbors should be zero. This leads to a disperse-less 
lower Hubbard band. As we see, Eq.(|97]) and (^8|) correspond to the standard model for the lanthanides if the upper 
transitions will not contribute to the cohesive energy and this situation arises when the energy of the upper transition 



A a2 3> £f- As known from electron spectroscopy, this is not the case even for rare earths EI However, as we have seen 
above, if the mixing interaction of the conduction bands with the upper transition is not equal to zero, the spectral 
weight at both, upper and lower energies, deviates from unity and, therefore, it is desirable to correct the localized 
model at least in this place. Let us introduce this correction. 

The exchange-correlation potential is a function of total charge density. Within the real-space representation of 



atomic-like wave functions this charge density is given by Eq.(lOO), where the Fermion GF (Taicfyw should be 



represented via the many-electron GFs, G. Within the approximation of dispersion-less lower transitions A ai we can 
write the structure of the XX-GFs as follows: 
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while all other GFs, G^ Xc \ C?( CC J and G^ cX \ do not have the lower pole due to the absence of hopping and mixing in 
this region of energy. Here a\ denote the lower transitions and 02 the upper ones. Inserting this formula into the 
expression for the charge density, Eq. (|100|) , we find in terms of retarded GFs 
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The last term represents the contribution of the core states into the spectrum of the Fermion-like excitations. This 
expression differs from the standard one used in ab initio calculations only by the numbers P a , which are present 
in the last term and all GFs. At zero temperature P a are determined mainly by the strength of the c-/-mixing, 
/-/-hopping and the proximity of A2 to the Fermi energy. In the last term the factor P a /[e' 3 ( Aii ~' i ) + 1] for the 
transition a = [7, T] gives the statistical weight with which the orbitals vq, v' contribute to the states I7), 

The expression for the total energy has the same ingredients, but the additional integral over coordinates gives the 
overlap matrixes and and matrix elements of the kinetic energy. Therefore, we find for the total energy 

1 f 1 

Etot -~2^J dc V(-A0 + l 

*{[Oj 1 L 1 ,j 2 L 2 uJ + t jlLl ,j 2 L 2 ]ImF h L 2tjlLx (uj + iS) 
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Here t is kinetic energy. 



IV. DISCUSSION AND CONCLUSIONS 



The ab initio calculations based on the LDA in DFT has become a quite reliable method for investigating the solid 
state. The materials with strongly correlated electrons, which forms a quite wide class (high-T c superconductors, doped 
Mott insulators, materials with colossal magneto resistance, heavy-Fermion systems and Kondo lattices and, different 
artificial structures based on small quantum dots) and are used in more and more applications, are traditionally difficult 
for ab initio calculations. It is well-known that these difficulties often are caused by the localized or quasi-localized 
electrons originating from open electronic shells of transition elements, lanthanides or actinides. Three methods found 
for overcoming these difficulties, LDA+U, self-interaction correction (SIC), and the complete localization model for 
lanthanides (where all /-electrons are forced to be core states) provide the same physics: they localize some of the 
electrons, making the potentials for conduction electrons and for localized ones, different. On one hand, since at least 
LDA+U and SIC use more than one potential, they go beyond the DFT scheme, where only one potential corresponds 
to a charge density in ground state. On the other hand, there always questions, How to avoid double-counting? What 
is "Hubbard 17"? Why is SIC applied to certain orbitals and not to other ones? etc., for which it is difficult to give an 
answer without considering the full problem (say, in the language of some of perturbation theories) and understanding 
what is included into LDA. An answer to first question has partly been found by Kotani li-3, who has shown that 
the LDA potential corresponds to a static random phase approximation. This simultaneously gives an answer to the 
question why the Hartree-Fock calculations give much less accurate results than the ones in LDA: LDA takes into 
account screening, while the Hartree-Fock method does not. However, this is only a part of the answer since, on 
the one hand, in modern calculations the exchange-correlation potential is taken from the (parameterized) results of 
Monte-Carlo simulations for homogeneous electron gas, which certainly, go beyond RPA. On the other hand, these 
statements and investigations are based on an analysis of either the homogeneous or weakly inhomogeneous electron 
gas. What can be done in the case of a strongly inhomogeneous electron liquid? The LDA+U and SIC methods 
are intended for this case and, as was mentioned above, they go beyond the DFT ideology and contains some open 
questions. As we see, the problem of understanding, in a formal language, what is included into LDA and what 
should be taken into account additionally in it to describe correlated systems are necessary steps in the development 
of the theory of solids. In the present work we made an attempt to obtain the physical picture analogous to the one 
generated by these three methods, but to perform such an analysis on the basis oLsplid grounds. Here we have used 
the perturbation theory from the atomic limit developed in a previous paper (refO) for the case of the generalized 
multi-orbital Hubbard- Anderson models and extended it for the case of inter-site Coulomb interaction. We have 
considered the easiest case, when a part of the electrons of the /-shell is well localized. The picture, generated by the 
perturbation theory is the following. 

First, under the word "correlations" we mean that the intra-atomic interactions are strong enough to make it 
necessary to calculate (fin) beyond the decoupled form (n) (n). Here, we use the exact relation (hn) = (n). This 
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immediately leads to the formation of a energy gap between the occupied and empty orbitals. Therefore, taking into 
account local correlations automatically leads to orbital polarization, the picture which is obtained in the methods 
LDA+U and SIC. The potentials for "upper" and "lower" orbitals differ from each other, roughly speaking, by the 
Hubbard U. 

Second, the language of the localized and dclocalizcd orbitals is not completely correct: since the orbitals are coupled 
by interactions to complexes, described by many-electron wave functions, the role of the Fermion excitations is played 
by the single-electron transitions between many-electron states. This is seen from the fact that switching on/off the 
mixing, (and/or hopping) between the upper transitions |r n ) — > |r„ + i) and the delocalized states, leads to a deviation 
of the upper and lower spectral weights from zero or unity and make these weights non-integer numbers. Then, the 
Fermion GFs acquire a many-pole structure, where each pole contributes with its own weight. In turn, these weights 
are determined by the population numbers (of the type Nr n ~ (^1^2. ..n n (l — n n -|_i)(l — n n -|_2)-..(l — h n))j where N 
is the number of orbitals in the shell) of the many-electron states involved. 

Third, the number of delocalized, correlated, orbitals (delocalized renormalized Fermions) does not obligatory 
coincide with the number of the single-electron transitions. 

Fourth, a part of the many-electron states which are composed mainly of the localized orbitals, often form a 
(quasi)degenerate subsystem. Redistribution of spectral weights between the transitions involving these states may 
be caused by small energies. 

Fifth, the spectral weights enter as factors to mixing, hopping, and Coulomb interaction, decreasing their strength 
and providing the correlation caused mechanism of band narrowing. Since these weights are determined by the 
strength of the interactions, the problem requires self-consistent solution which may reflect the competition between 
the Hund-rule interactions and the cohesive energy. 

Sixth, the closed equations for the GFs of the Hubbard operators in terms of functional derivatives contain two 
types of contributions. One type is responsible for the description of effects from kinematic interactions originating 
frora-piixing and hopping ( an example of the physics involved, the correlation mechanism of derealization, see in 
Ref.EJ), another type generates an analogue of the WCPT for the Fermion GF. The series for the many-electron GFs 
obtained by neglecting the contributions from the kinematic interactions is isomorphic to the series for the single- 
electron Fermion GF in WCPT. Using this isomorphism, and comparing the Migdal-Galitskii formula for total energy 
of the non-correlated and correlated systems as well as the Sham equation for the exchange-correlation potential, we 
came to the conclusion that the form for the exchange-correlation potential used within weak-coupling theories can 
be used in the theory of the correlated systems also, if one will use the mixing, hopping, the inverse of the overlap 
matrix and the expression for charge density renormalized by many-electron spectral weights. The recipe for concrete 
calculations and example of such calculation is given in future works. 

Seventh. The model Hubbard U can be expressed in terms of energies of transitions between many-electron states. 
Consideration of screening of these transitions by Coulomb interactions shows that the lower transitions are renormal- 
ized much weaker than the upper ones. This makes calculations within one parameter, U mod , somewhat confusing, 
since they do not reflect the true picture. 

Let us discuss this question in a more details, and make an attempt to elucidate what actually should be understood 
under Hubbard U from the point of view of this many-electron theory. A commonly used definition of the Hubbard 
U has risen from the following speculation. The energy of the n-orbital atom modeled with the Hubbard model is 

E n = ne d + ^n(n - l)U mod . (102) 

From this one finds that the energies of single-electron excitations are 

A 1 =E n -E n _ 1 =s d + U mod (n-l), (103) 

A 2 = E n+1 - E n = e d + U mod n (104) 

and, therefore, 

U rnod = A 2 - Al (105) 

The parameter £0 can be taken from the requirement that the effective attractive potential is able to hold n electrons 
(this is determined by the valence of the atom) and, therefore, £q = — (n — l/2)U mod and, in turn, Ai = —U mod /2 
and A2 = U/2. Thus, Ai and A2 play the role of the centers of the Hubbard sub bands. Hubbard used the Coulomb 
matrix element of density-density type for jj mod _ According to the common opinion, it has unreasonably large values, 
between 20 and 30 eV for a d-shell, and, therefore, is considered irrelevant. In LDA+U calculations "reasonable" 
values of jj mod for transition metals are near 5-8 eVH, while in recent calculations, based on dynamical mean field 
theory, in order to obtain a reasonable agreement with the spectroscopical data for iron the authors had to take a value 
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of about 3 eVH3. There are, however, no suggestions what mechanism that is responsible for such a drastical decrease 
of U from 20 to 6 (or even to 3)eV. The common believe is that derealization plus screening, particularly, due to 
interaction with some other bands, may provide such a decrease. It is easy to see, that at least, the density-density 
interaction with other bands does not change the picture. Indeed, let us imagine, that besides, say, a d-shell the atom 
also contains another shell with m atoms. If we add the term U c d,mn to the energy, it changes Ai and A2 

Ai Ai + U cd m, A 2 => A 2 + U cd m, (106) 

but does not change U mod : A 2 - Ai = U mod . EarlierB we also investigated the influence of the derealization on 
jjmod y j a ^normalization by ( kinematic) mixing interaction, it gives a decrease of at most 10% under favorable 
conditions. Therefore, either other mechanisms should be involved or the description in terms of one parameter U mod 
is oversimplified. From our consideration given above we can conclude that both assumptions are relevant. 

First, it is important, which experiment we want to compare with. It is well-known that for ab initio DFT LDA 
calculations the excited states with large energy of excitations are difficult since they usually involve non-single- 
electron excitations and the latter are not included into the effective single-electron Kohn-Sham scheme (even with 
different corrections). However, in low-energy experiments the lower Hubbard band is not involved. Therefore, for 
the description of such experiments one is interested only in the properties of the upper Hubbard band. The model 
calculation, which are trying to work with energies separated by a large scale, usually does not take into account self- 
consistent changes of the charge density. This involves the strongest interactions existing in the system. Therefore, 
it looks like only methods which combine the ab initio methods with the many-body theory might be able to treat 
such problems. Within our many-electron approach the important parameters of the system are Ai and A 2 . They 
are connected with the model parameter Hubbard U , however, they have, obviously, many-electron nature. 

Second, there are two reasons why it is very important, where the levels Ai and A 2 are situated with respect to 
all other bands and the chemical potential: a) this influences the distribution of spectral weights; b) the efficiency of 
screening depends on it. Actually, the lower the transition is under the Fermi level, the less the Coulomb interaction, 
via which it interacts with other electrons, is screened (see the speculations about the one loop contribution to the 
polarization operator, Eq.([78|)). 

Third, the way of calculation, obviously, depends on the reference point which interactions are included into 
A*f\ If we start with the LDA picture with a certain number of localized electrons, we have self-consistent 
/-wave functions, from them we can construct of the wave functions of the many-electron states, and calculate the 
energies i , E^\ E^ +i without any relaxation of charge. The differences of these energies provides A^\ . It 
should be emphasized that in this approach the LDA based calculation is used only as a generator of a convenient 
minimal set of wave functions for calculating E^\ 1 , E$?\ . The energies themselves should be calculated on 

the basis of the many-electron states, Then, one can switch on mixing and hopping and take into account the 
screening and the renormalization of A^, A^ "* via the self-energies of the many-electron GFs. Let us start with 
the lower transition. If we denote the Coulomb interaction between the localized and delocalized electrons as Uu, 
delocalized-delocalized as Udd, it is clear that the largest contribution to self-energy for the renormalization of Ai comes 
from the graphs which contain the lowest possible number of interactions Um, because, as we have seen in previous 
section, each loop constructed on Uu involves the loop with large denominator of the order of A 2 — Ai = U mod . 
Therefore, the largest contribution to the self-energy with density-density matrix elements Uid, renormalizing the 
lower transition, contains Um only twice (see the graph in Fig. [l4| and Fig. |l5|). Strictly speaking, the equations for 
the lower (and for upper as well) GFs always form a system of equations, 

a i a l' (i, ,\ _l_ /^ a l a 2 (:, .\ rvvff 



G ai3 i' (iw) = Go iai ' (tw) + Go lQ2 M{£ 



a 2 a 1 



a 2 69a 5 1197,68 (iv - iwi) 



x^ 3 78a 4 G a5a3 ^i)]}G a 4 a i'(zc;), (107) 

where £^ Q4 is the Hartree contribution and n is the full polarization operator. However, if we neglect the differences 
between the levels Ai, we can write the equation for renormalized Ai: 

Ai(iw) = A? + aT }^ t%i69pi 1^97,68 - i^i) 

0)1 

XW ai78 a 1 G aiai (^i), (108) 



where a takes into account the number of internal transitions, 05 in Eq.(107). The situation with the upper transition 



is more complex since in this region the bands are formed. For this reason it can be defined via the static dielectric 
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per meab ility as a solution of the equation with Hartree and, screened by static exchange interaction, i.e in the 
Eq.( |l07|) using instead of the second term with £ _1 (0) and F and v replaced by G and v. As we see, the self-consistent 
magnitudes A2, Ai become frequency-dependent already in the first steps of their dressing by electron-hole excitations. 
This makes such calculations quite cumbersome and practically difficult. The other possible way would be to calculate, 
as described above, the non-self-consistent values of A 2 °\ A^ and to treat them as the input parameters for further 
LDA calculations, where they will be renormalized as centers of bands usually do. 

In conclusion, we have reported the following results. The perturbation theory for correlated electron systems 
developed in previous paper here is extended to the case of Coulomb interaction. It is shown in the level of closed 
functional equations for the many-electron Green functions that all the contributions in SCPT can be separated 
to two classes. One of them forms the approximation of renormalized fermions, ARF, does not include kinematic 
interactions, and is fully analogous to the WCPT for single-electron GFs. Another one includes them and, therefore, 
contains specific for CES effects. The ARF is used for formulation of simple recipe for extension of the LDA to DFT 
to the case of strongly correlated electrons. The application to rare earths is given in next paper. 
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APPENDIX A: THE COULOMB TERMS IN EQUATIONS OF MOTION FOR GREEN FUNCTIONS 



Class 0. This class belongs to the ones considered above. Nevertheless, for completeness we give the detailed 
formulas. The correction to the EM for the GF {Tci(t)rf{t')) u (77 is either a c or a X operator) from the interaction 
between the conduction electrons is 

xCTc^t (;'))„], (Al) 
while the GFs generated by the commutator [X, Vo] can be expressed as follows 

(TZ^(t)cl(t)c 4 (t)c 5 (t) V Ht>)) u 

= lim [{(TZf^u+i *. . ) 

x(Tct(t)c 4 (t)c 5 (tW(t'))u], (A2) 
= lim [((T^ 1 (*!))„+«— |— ) 

x lim , „«T4(*2)c 4 (t a ))u + *: 



x(Tc 5 (t)rf (*')>„]. (A3) 
Below, for brevity, we shall denote the time limits by upper index " plus" : 

lim R(t 2 ) = R(t+). (A4) 

Further, 

{Txf i (t)c i (t)c 5 (t)r 1 Ht')) u 

= {{ TXf 1 (t+)c i (t+)) u + i —J — ) 

x(Tc 5 (trf(t>)) u . (A5) 
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Then in the Hartree approximation we have to neglect the contribution from this derivative. We shall use the notation 
( j(o)£<( c 'j) f or j-jjg right-hand side contribution to the EM for corresponding GF, not for the GF itself. Thus, for the 
GF (l/i)(Ta(t)rf (t')) u — we have the standard expression 

s( o) G (*n) = V^^TcUt^cSiViTcSWit'Vu, (A6) 

where Vi c 3 c 4 c 5 c = Oi^ 2c v 2 c 3 c 4, c 5 c - The index (0) means that the correction comes from the class 0. This is, obviously, 
the standard Hartree correction. The exchange-correction comes from the derivative. Next, 

x(T^{t + )c A (t + ))u(Tc,{tW{t>)) u 

x{rxl{t + ) Ci {t + )) u {rcSW{t'))u. (A7) 

Class 1. The first GF arising from the c-operator is (TZc^ccr/^} which is considered above. The other functions 
coming from [c\, V{\ are 

(TX%(t)c 4 (t)c 5 (tW(t'))u 

x(Tc 5 (t)rf(t')} u , (A8) 

and 

{Tc\{t)x«;{t)cSW{t')) u 

x(Tc 5 (t)rf(t')) u . (A9) 
The functions which come from [c±, V\] we rewrite as 



s 



x[<T4(t+) C4 (t+)) u + »_ Fy ]<Tc 5 (t)^(t')) u ] 



" W [2 c 3 c ]4 / 5 c (^2 t ) a2 °l c 3 c 



34^+V 

6 



x[(TX%(t+)c 4 (t+)) u +i " }(Tc 5 (t)rf(t')) u 

+' y [2 c 3 c 4 / 5c]^4 40 l c 1 2 c 

x[(Tct(t+)X%(t+))+i 6 ](Tc 5 (t)r)Ht>)) u . 

0U 3c,i4.aA l I 



Therefore, the Hartree contribution is 

5 = » [U ]^ t ^) i2 °4/^r i! ( T 4(' + )>« 

-« [ 2 c 3 c ]4 / 5 c (/2 t ) S2 Or c 1 3 c [(^ 2 (t + )c 4 (i + )} U ](Tc 5 (i)^ t (i / ))« 

+V3 c 4 / 5 c ]/4 Q4 or c 1 2 c (r4(t+)x;;(t+))(r C5 (t)r,t(t'))„. 



(A10) 



(AH) 
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The contribution to the GF (TXrf) is 

- «[ 2/3c ]4 c 5 c 4 ia2 (/2 t ) a2 ^ 1 ,J(^4( <+ ))^ 



5 



x [<T4(t+)c 4 (t+)) u + i gu J {t+) ](Tc 5 (t)vHt'))u} 

+^ / 3 c ]4 c 5 c Or / 1 3^r a3 (/ 2 t ) a2 (/3 t ) a3 4, 2 [(^4( i+ ))« 

(5 1 

x [(TXff (t + )c 4 (*+))„ + i 6 ] (Tc 5 (t)qt {t ')) u 

+ U [2 / 3 C ]4 C 5 C ( ^1/ 1 3 C £ C 1 3 (/2) a2 (/3) a3 ^'ii2' C b 2 ^ 

x[(TXi(t+)c 4 (t+)) u + i — J- ]{TcSW{t')) u 

+^3 c[ 4 /5c] o- 1 L, 2 jr(4) a2 e r a2 K^4(* + ))«+ i ^^] 

x[(T4(t+)X%(t+)) u +i 6 }(Tc 5 (tW(t')) u 

0U 3 c j4.a,4,\ > 

+^ c 3 [4 / 5 ] 7 1 l3,3 c /4°n4) S3£ r S3 {^, 4 [<^4(t + )^ 4 (* + ))« 

x «t4 (*+)>„ + » ^/ (t+) ]<^(*)'? t (* , ))t. 

+^/c^[<T < j(t+)^ 1 (t+)) tt + i — ^ ] 

x(Tc 5 (t)r,t(t'))4- (A12) 
The corresponding Hartree part is 

= ^ /3c] 4 C 5/r S2 (/2) S2 ^^( T 4(* + ))« 

x<T4(t+) C4 (t+)) u <T CB (t)^t(f)) u 

+«[ 2/ 3 c ]4 c 5 c r / 3c e r a3 (/2 t ) a2 (/3 t ) a3 ^, 2 (^4( <+ ))« 

x(r^(t + )c 4 (t+)) tl (T C5 (t) ?? t(t')) tt 

—1—7) O^ 1 p a l a 3/ J-t^a 2( ' ft\«3X . . K °-2i 

^ U [2 f 3 c ]4 c 5 c lf^ct V-/2J U3J "3132^1 

x{TXl(t + )c 4 (t+)) u {Tc 5 (tW(t')) u 

+«2 C 3 c[ 4 /5c] Oi;L, 2 j4 a4 (4 2 ) a2 4 ia2 ( r 4(* + ))« 

x<T4(t + )^(t+)) u <T CB (t)^t(f)) u 



+v n Q r.. Kc] 07 1 i M8l3e /r(/X) B3 er B H^^4(* + )^(t + )>t. 

x(TZ) i (i+))„(Tc 5 (t)r ? t(t'))„ 



+<5 J . 1J . 4 / S ^<T4(t+)^ 1 (t+)) t ,<Tc 6 (t)^(f)) u }. (A13) 
CZass 2. Here the corrections to the equations of motion are 
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s_ 

+hh /r (fir 2 iftro^T 2 \(tz% (t+)) u + 1^^+ 



*[(TX$ 2 (t + )c 4 (t + )) u + i u+ , ]{TcSW(t>)) u 



x[<TX°'(t+) C4 (t+)) u + » / (f+ A (Tc B (t) V H^ 
x[<TXff (t+)c 4 (t+)) u + » / r . + J <T C6 (t)»?t(0) u } 

0U 72a2,4 c V r J 



xl(TX%(t+)c 4 (t+)) u + i 5 }(Tc 5 (tW(t% 



(A14) 



+<w (fir 2 (flf^r 2 [(tzi (t + )) u + i 5uZ ^ t+) : 

x[(T4'(t+)c 4 (t+)) B +i ^x/ 4 (f+) ](^ 5 W^(*0) u 



x[(TXff(t+)c 4 ft+)) M + i / ](Tc 5 (tW{t')) u }. 

0U 3ia 2 A a \ l > 

The remaining contributions, in this class, comes from the processes of transforming two c- into /-electrons 

s ^G(f = \v Wc ^hMfiro^ s n^r 

x[(TXl(t+)c 4 (t+)) + i * {t+) ]{TcSW(t'))u 

+s j2j3 f?(flr 2 (flr 3 o^r 2 

x[(4( i+ )) +*^pry][(r^(t+) C4 (t+)) 

-s j2j jr(flr 2 (flr 3 o^r 3 

x[(4( i+ )) +*^pT)H(r^(t+) C4 (t+)) 



(A15) 



+ ^ TT+i KTcsitWV))*}- (A16) 



and 



S^Gj^ — 2 W 2 / 3 / 4 c 5 c {^j2ji^2i3 {flftr 5 ^ 71 

x[(TXl(t+)c 4 (t+)) +i 6 ](Tc 5 (tW(t')) u 
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+ ^2jl ^32 33 (fl)"" 2 (fl lCt2 

-^ 1 4, 3 (/ 2 t ) a2 (/ 3 t ) a3 4 ia3 

x[( ^ (<+)> +l ^7FT H(r ^ 2(t+)c4(t+)> 

(f+ , WcSW{t')) u }. (A17) 

Thus, these two corrections differ only by factors of the overlap matrixes coming from the commutator of c-operator 
with interaction. Class 2, 3e). Indeed, a c-operator, located on site j\ "feels" strong correlations on other site j'2 only 
due to non-orthogonality, i.e. due to the /-component which it contains on the site ji. Remaining terms in this class 
come from the interaction 3f, where two /-electrons are transformed to c-electrons: 

x{Tc\zl{t)r,\t')) 

+6 jd5 [(TcUt+)Xt(t + ))+i J ] 

0U 3 c ,j 4 a4.\ l I 

x(TX^(t)rf(t'))}, (A18) 



and 



x(T(clzl(tW(t')) 

+hh [<4 (* + )) + i 5U z^ t+) \ K^4(* + )*r (t + )) 

+i 5lJ J m ]{TX?{tW{t>)). (A19) 

We see that new types of GFs, (Tc^Z^^t)^ \t')) , have appeared which cannot be decoupled in a standard way. 
Expectation values of the type (Tc^it)^ (t 1 )) or (TZ^it)) are equal to zero in the equilibrium non-superconducting 
state. In the external fields J dtU v (t) Z v (t) , they do exist. Therefore, we should keep these expectation values while 
the fields are not equal to zero and after calculation of the derivatives S(T Z] 4 (t)) / '8U 71 '(f) and 5(Tc\ (t)c\t')) /5U cc (t') 

, 5(Tcl(t)X a (t')) /5U Xc (t'), etc., We can put them equal to zero. Usually these graphs give small contributions 
since they involve the gap ~ U in the denominator. However, these contributions describe the an effective anti- 
ferromagnetic interaction between (quasi)localized electrons and, therefore, without their contributions it is difficult 
to obtain localized antiferromagnetism at all. Thus, we can use, for example, the decoupling 

(Tc\zl{t)n\t')) = [(Zl(t-)) +» su J — }(Tct(t)^(t')). (A20) 

The direct derivative of the zero inverted GF, (Tc^t)^ (t'))^ 1 , is equal to zero. However, it is always connected in the 
equations of motion with the operators cZ r > that gives the non-zero loop correction from (Tc 3 (t)c\ (t'))(TZj 4 (t')Z^, (t)) . 

Class 3. The terms, containing three /-operators may have the operators, belonging to the same site, or to different 
ones. If they belong to the same site, we can write them as follows 

^2/3/ [4/5 c ] ^hh &h33 (flflfiT 2 Xj 2c 5 

+^2 c s f] 4 f5f ^hjJuh {fthhYcUl ■ (A21) 
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They, obviously, also belong to mixing interaction and have been included to it in Ref.Q. 

Next are the terms containing two /-operators, one of them belonging to one site, and another one to other site. 
They can be written as follows 



x{o^ 2f (f 2 retr k4 2 ( t+ » +* su z S (f+) 1 



' ' Z t( t+ ))+ ^ T J ,, + s ](Tc 5 (t)r,Ht)) l 



Then, for the GF (TX 81 (t)^ (*')) we have 



0U j2a 2 J 3 b\ l ) 
+ U [ 2c 3 / ]4 / 5/^34,3(/3 t /4) C (/5) a5 {0r c 1 2 J(4( <+ )) 

-or„ 1 3 / (/3) a6 ^ a [<^(* + )^i,(* + )) 

V^ 1 ™ ,(0,} ' (A22) 



><{^4f 2 K^(* + )) + ^^/( F )] 



x[(4 3 ft + ))+ v t/ / (t+) ]< Tc 5 (^ t ft / )) 1 



0U 32a 2 ,jiby T ) 



+ V [2 c 3 1 }i f b 1 ^3i33^3s33(flf'i)^{fb) a& 

x {or, 1 ,. (4re°r [(4- (<+)) + i * ] 
><[<4(* + ))+ i ^|^K TX "(^ t (*'))« 

-«i li ,er ,fi [^4(* + )^ 1 (* + )> 



Ctass ^. Since the terms, which contain four /-operators and belong to the same site, are included into the definition 
of the on-site energies E-p of the terms T, we start with the density-density interaction between different sites 

§(4)Qff — -v 2f3f4f5f dj 5 j 2 Sj i j 3 Sj 2 j 3 (flf 5 )^ 2 (flf 4 )^ 3 

x{or c 1 2/ (/ 2 ) Qe c ?2 [(4( i+ )) 

+l «7^(F) ](r ^ (<)?7t(< ' )>tl (A24) 
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+i su J {t+) ](TXl(tW(t'))uh (A25) 



and 

(t+) }{TXl(t)rf(t>)} u (A26) 

+5 3l jA li3 \{zf 2 {t + )) 

^ m z\ t+) W^ 3 {tW(t')) u }. (A27) 
At last, the corrections to the EM from single-particle potential are considered in details in Ref.0. 



APPENDIX B: SPECTRAL WEIGHT TRANSFER: AN EXAMPLE OF TWO LOCALIZED ELECTRONS 

IN THREE-ORBITAL ATOMS. 

The LDA+U and self-interaction correction methods are used in DFT for creating an orbital polarization, which 
allows to describe localized and delocalized electrons with different potentials. In this section we illustrate the 
mechanism of orbital polarization generated by strong intra-atomic interactions considering the example of an ion 
with 2 localized f-electrons. In order to make consideration more transparent we consider a hypothetical 3-orbital 
atom with spin less electrons. This electron shell can be described by the Hamiltonian: 

H d = H° d + H u dd = £on m + U h rnh m ' , (Bl) 

where m,m' = 1,2,3 and h m = dl n d m . The parameter Eq describes the attraction to the nuclear potential which 
provides the ground state with two bonded electrons 

= => e = -hj. (B2) 

n=2 1 

The Hamiltonian has the following eigenstates and eigenvalues 



1 / \ dE n 
E n = e Q n + -Un(n - 1), — — 

2 dn 



n 


state 







|0) = |000) 





1 


4|o) = |100) 


2 U 


1 


4|0) = |010) 


-It/ 


1 


4|o) = |001) 


2 


2 


44|o) = |no) 


-2U 


2 


44l°) = l 011 ) 


-2U 


2 


44 1°) = koi) 


-2U 


3 


444\o) = |in) 


-^U 



Introducing the Hubbard operators 



■■ {\ — 5ij)X lJ +8ijh % in terms of these states we diagonalize the Hamiltonian 
- 2U[h lw + h 011 + h W1 }. (B3) 
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The Fermion operators in terms of Hubbard ones are 

di = (OOOldillOO)^ 000 ' 10 ^ (010K|110)X 010 < 110 
+ (001|d 1 |101)X 001 < 101 + (011|di|lll)X 011 ' m 

_ j^000,100 _|_ j^-010,110 _|_ jj^-001,101 _|_ ^-011,111 

d 2 = X Qoa < 010 - X W0 - 1W + X 001 ' 011 - X 101 - 111 , (B5) 

^ _ ^000,001 _ ^-100,101 _ ^-010,011 _|_ ^-110,111 ^gg^ 

Correspondently, the Fermion GFs are 

F n (iw) = <Td 1 (T)4(T , )) ia) 

= -(T(X om ' wa + x 010,110 + x 001:101 + x 011,111 ) 
* (x ioo,ooo + x iio,oio + x ioi,ooi + x iii,oii )t/) ^ 



I(TX 000 ' 100 (r)X 100 ' 000 (r'))^ 

+ \ { T(X 010 ' no (r)X no ' ow (T% u 
+ 1 (tx ooi,ioi (t)x ioi,ooi (t / )) . w 

+ 1 (tx oii,iii (t)x iii,oh (t / )) . w 

. N 00 o + N 10 o . Now + N 110 

— l = l- 



. Nqqi + Nioi . Non + Nm 
-i -, i- 



1/ -i V -I 

iu) + o(7 + fi iu — yU + \i 



(B7) 



The non-diagonal GF like \(T(X 0Jl (T)X IIlj2 (T')} iuJ = due to the diagonal form of the Hamiltonian in this 
representation. In the same fashion 

F 22 (iw) = (Td 2 (T)4(r')) iu 

. iVooo + N w . iVioo + Nno 
i- 



- N n _ iVioi + i 
iu + \U + /i iuu — }-U + ^ 



.AW+Noii .N W i+Nm 
~ % - — lrr , % - TT1~, — ' ( B8 ) 



F 33 (iu) = (Td 3 (T)4(r')) iu 



. N oo + N 00 i . ^Vioo + Niqi 

-l 5 l- 



iu+^U + n iuj + \\J + fi 

. Now + Non . Nno + Nm 

— i -, i t— . (B9) 

iuj + j^U + fi iuj — T^U + ji 

The population numbers in the numerators for the non-interacting atoms are just the Gibbs weights, 

iV oo = 1/Z , iV 100 = No W = Now = e^+^/Z , 
N 1W = Non = N W i = e^ +2U ^/Z , 

Zo = l + 3e^+i u 1 + 3e /3(2p+2C/) + e/9(3M+§iO (B10 ) 

The chemical potential has to provide that the total number of electrons in the ion equals two. The equation for the 
number of electrons in this ion, 
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n d = Y,(4d t ) = 2, 

i=l 

or written in terms of the many- electron population numbers, 

n d = • W 00 o + 1 • (A^ioo + W i + Who) 
+2 • (N 1W + JVon + W 10 i) + 3 • Win 

3 e /3(M-|tO + 6e /3(2 M +2(7) + 3e /3(3 M +§tO 



(Bll) 



1 + Se^+f* 7 ) + 36/3(2^+2(7) + e /3(3 M +|C/) 



(B12) 



can be fulfilled for \i in the interval: Ai = E 2 — Ex 
population numbers fulfill the condition 



■¥ < A 1 < t = A 2 = E 3 — E 2 . In this interval of \x the 



Wrao = Win = Wioo = Woio = Wooi = 0. 
Since the sum of all population numbers is equal to unity, 

l = ({L,f!})=J2 N r> 



(B13) 
(B14) 



and, also, 



all three zero Fermion GFs, 



Won = Wioi = Wi 



10 



Wi 



10 



+ 



W 



101 



Won 



ioj + ill + n ioj + ill + n iu> — ill + /J, 



(B15) 



iF° 2 (ilu) 



Woi 



Wic 



(B16) 



W 



101 



IUJ+W + \i 



+ 



Won 



+ 



Wn 



become equal to each other within the statistical description, 

^Fl lM - iF&M - iF° 3 {iw) 

2/3 1/3 



iuj + W + fi 



+ ■ 



(B17) 



(B18) 



As seen, when the chemical potential is within this interval of energy, 2/3 of the spectral weight of each orbital is 
concentrated at the energy ui = —\U — /i, while the rest of this weight is placed, by the sum rule, at u> = — fi. 

From the Fermion GFs we find that the total spectral weight in the lower and upper poles is the same as within 
the single-electron picture, 

^Fl lM + iF° 2 (iuj) + iF^iLu) 

2/3 1/3 



+ 



IUJ + \U + jJL IUJ - \U + jJL 



ioj + \U + [i iui — -^U + fi 
and the same number of particles for the chemical potentials in this interval of energy 



(B19) 



n d = nx + n 2 + n 3 
= 3 2/3 



+ 



1/3 



e /3[-2U-A + I e P[-3U/2-n] _|_ 1 



= 2. 



(B20) 
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This explains why the phcnomenological recipe of the DFT LDA calculations with some fixed number of the localized 
electrons (in the case under consideration, two) works well: at large values of U the contribution of the upper pole is 
not taken into account and when the states are degenerate, there is no difference, which of the orbitals of the non-filled 
shell that are occupied. However, it is instructive also to go further and consider the case when the degeneracy of 
these states is lifted away. 

Let us imagine that some interaction pushes the energy of one of these orbitals below the other two (e.g. Hund 
rules, or long-range order). Then, we have to shift down the energy of, say, the first orbital, — > £o — h e ff, arid, 

(2) (3) 

£q — > £0) e o ~ * £ o + heff m tnc Hamiltonian. The population numbers of the levels without electrons and with 
three electrons are still equal to zero because h e ff <C U. Although the energies of the three orbitals which share two 
electrons differ from each other by the magnitude h e ff <C U, at low enough temperatures <C h e ff the population 
numbers become not equal to each other, 



f - 
1 



Nno i + ePK,, + e -Ph eff ~ 



N wi 1 + ef3heff + e -0h eff ~ °' 
N ou = 7TT m, — ~ 0, (B21) 



this leads immediately to a redistribution of the spectral weight 

• no / ■ \ Niio A^ioi 

+ 



iu: + ±U + [i iiv + ±U + h e ff + [I 
Nou 



iu>- \XJ + h eff + [1 

(B22) 



2 l 

1 



iF° 2 (iu) 



Ano , -Aon , Aioi 

(B23) 



iuj + \XJ + fj, iuj + \XJ + iuj - \XJ + fi 
1 



iw + \U + n 



Aim Ar 



iu + \U + [i iuj + \U - h e ff + n 

A 

+- 



2 l 

Ano 



iuj - \U ■ 



iu - \U + fi' 



(B24) 



Thus, in this case, the spectral weight of the two localized orbitals is concentrated near the lower poles, iuo = 
£(110) - £(100) ~ -\U, while the spectral weight of the third orbital is concentrated at the much higher energy, 
iu ~ +\U. This phenomenon is the orbital polarization, and it is caused here by the combined effect of the effective 
field and the sum rules. As we see, actually, the phenomenon of the orbital polarization is caused by two (in our 
case) different interactions which involve very different energy scales. First, the strong Coulomb intra-atomic repulsion 
compensates for the strong attraction of electrons by the nuclear potential to the extent, corresponding to the valence, 
displayed by the particular ion in the particular surrounding. The latter, should be taken into account via a Coulomb 
renormalization of the energies of the single-electron transitions and it is considered in the discussion and conclusion. 

Second, a much smaller energy than the "Hubbard U" and Hund- rule integrals, like an energy of magnetic ordering, 
or, crystalline electric field, can provide a further polarization of the electronic shell. The Fermion GFs have purely 
single-electron form with location of orbitals at very different energies and, as we have seen, it is the intra-atomic 
interaction plus the sum rules for the many-electron spectral weights that provide their effectively single-electron 
form. Up to now we discussed only fully localized states. 
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Let us now discuss the following questions: 

1. How does the delocalization of the non- filled orbitals arise and how does the the spectral weights change? 

2. What approximations should be done in order to be able to separate the calculations into core and valent ones 
as it is usually done in ab initio calculations? 

To consider these questions, we have to add to the Hamiltonian the mixing interaction with some conduction band 
and the hopping term: 

H = (Hd + Ti^d) + T^c + 'Hmix + 7~ihop', 

K = E £ 4c fe ; n hop = J2 ti^Ujui 

k (ij)v 
k 

Here we have chosen a simple form of the hopping, i.e. hopping is allowed from the orbital v in one site to the same 
orbital in the neighboring ion, in a real lattice it is often not so. 

The first immediate conclusion comes from the form of the Fermion GFs, F. Only a part of the /(d)-electron shell 
is delocalized since, at the energies uj ~ —\TJ, the hopping matrix elements tn, t\2, ^13, £23, ^22 are exponentially small 
(ipf ~ exp (— J dxy/2m[U(x) — Ai]), where U(x) is the potential separating the /-electrons belonging to neighboring 
atoms) and the transitions Ai = E2 — Ex = —\U remain localized, while the A2 = £3 — £2 = may be located 

not very much higher the chemical potential and the magnitude of the hopping matrix elements £33 may be noticeable; 
the mixing interaction is effective also only with the orbital which energy crosses a band. Then it looks like if we are 
back to the standard band picture, but only for the "upper" orbitals; therefore, the expected (index "exp") equation 
for the Fermion GF for the third orbital is 

F-Jp-M^) = I^sMr 1 - t 33 (k) - ] ] k ] }\ (B25) 

and the important physical consequence expected from this—form is that the tails of this upper f(d) band may go 
under the Fermi level and contribute to the cohesive energyEx Below we will see that this single-electron approach 
misses the important transfer of the spectral weight. 

Let us consider the problem a bit more accurately. First we find the GFs for the transitions involved, the equations 
for the population numbers and, make conclusions about the spectral weights and the approximation needed to be 
done in order to make the problem as close as possible to the single-electron picture. First, since only the term 
^33(^)^3^3 is effective, the main role in the formation of the collective state will be played by those transitions for 
which the matrix elements d^ jr ^ = (7^3 |T) is non-zero. It is convenient to introduce a short notation for them 



2 = [110,111], 2 = 
1 = [100,101], 1 = 
I' = [010,011], F = 

Then, the system of equations for the GFs within the Hubbard-I approximation. 

T) H 1 IA G HIA = P, (B27) 



(which is sufficient for our target here) looks as follows 



[111,110]; 
[101,100]; 

[011,010]. (B26) 



n 2 - A 2 A 2 A 2 \ / G 21 G 21 ' \ 

Ai n 1 -X 1 -Ax G 12 G 11 G 11 ' 

Xv -X v Qy-Xv J yc 1 ' 2 G 1 ' 2 G 1 ' 1 ' J 

/ P 22 \ 

= P 11 . (B28) 
V P l ' v J 

Here we we have used the notations for transitions 2, 1, 1' used above (don't confuse with the orbital numbers): 
G 22 = -i(TX 2 X 2 ) lLU = -i{TXl U0 > ni lXl in > n % u , etc. 
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flj = iuj — Aj; Xj = X(k, iui) 



|V(A;)| 2 



iuj - Sk 



Pj, J -2, 1,1', 



(B29) 



Then, 



G 



HI A 



\D H \ A \ 



I', I'": P 2 =N U1 +N W1 ; 
Pi = N wo + N W1 ~ N W1 ; 

Pv = N 0W + N ml ~ N ml . 



SivPi -X 2 fl v P 1 -A 2 f2iPi- 
-Aifivft $2i'Pi Aif2 2 Pi' 
-A1/O1P2 -Ai/n 2 Pi $2iPi' 



(B30) 



(B31) 



where 



l-D 



— Q^Qj fl^Aj AjOj. 



(B32) 



Let us write down the Fermion GFs via the Hubbard GFs. We start with the Fermion GF of the dclocalizcd "orbital" 



Fyy A = d v {d\,) b G ab = G 



f \b r^ab 



G 



21 



G 



21' 



-G 



12 



G 11 + G 



11' 



G 



l'2 



- G 1 ' 1 + G 



p 2 n 1 n v + n 2 p 1 n 1 



i'i' 



Q 2 £llQy — X 2 ^li^ly — ^ 1 ^ 1' — ^2^lAl ; 



(B33) 



and inspect the form of this GF in two regions of energy, iuo ~ A2 and ioj <~ Aj. In order to see a rough structure, 
let use the fact that the magnitude of the effective field, h e ff, is invisible on the scale of energies ~ U, therefore, we 
can put f^i = f2i>; and 

P2 + P1+ Pv = {Nuo + N in ) + JV101 + JV011 = 1, (B34) 
(with the accuracy N W o = N i = ^001 — 0); then, we see that the approximate GF for the third orbital is: 



pHIA _ 



P£h + (1 - P 2 )»2 

o 2 Oi - A(i - p 2 )n 2 - AP 2 fii ' 



(B35) 



Then, we have a third-order equation for the spectrum. In order to make the further formulas readable we switch off 
the mixing since this does not change the principal features of the model. From the poles of this GF we find the band 
spectrum: 



E 



1 



k - ^(Afc + Vk); v k = v /(C/ + A fe ) 2 -4A fc [/(l-P 1 ). 



(B36) 



This formally coincides with the spectrum of the s-band Hubbard model in the Hubbard-I approximation if we take 
into account that, in the form of the spectrum, we used the choice £0 = — 3U/2 explicitly The GF of the delocalizcd 
orbital can be written then in the form: 
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El 



+ 



where 



X k - U(l - 2P 2 ) 



vl = ± 



(B37) 
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Now we have a look at F22 and F33. Since t 2v = V 2v = 0, and the matrix elements d\ in the expansion 



7 HIA 

11 



dl{d\) b G 



t \bs~iab 



(B39) 
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are equal to zero for the transitions a = 2 = [110, 111]; 1 = [100, 101]; 1' = [010, Oil]. Therefore, 



T?HIA _ 

N 1W + N W1 Nou + Nui 

22 iu + U/2 + h e ff iuj-U/2 + heff' 1 ; 

Nn + N ii , N wl + Nm 
F " - iu + U/2 + iL u-U/2 ' (B41) 

i.e. the only difference with the case when the hopping and mixing are absent is that the population number Nm ^ 0. 
Indeed, we find from from the GF G 11 : 

Nm = (Nm + Ni W )J2i4.f(Ek - A*) 

k 

+b 2 k f(E+ - »)] (B42) 
= (JVin + JVno)e ? 0. (B43) 



Here, 

a k = 



2 1 



2 



1 U + \ k - 2\ k {\ - P 2 ) 

Vk 



,b 2 k = l-a 2 k . (B44) 



We consider the situation when the bottom of the conduction band is above the energy of the lower transitions, 
^bottom > Aj,Ap, an d the energy of the upper transition is above the Fermi energy, A 2 > £f- Since the level 
A 2 = U/2 before switching hopping and mixing was empty, the band being transformed into E^ changing is filled 
only slightly and, therefore, 9 « 1. Within the accuracy we need this makes it sufficient to replace the sum Q(P) 
by 0(p(°)), i.e. to find it without self-consistency, using non-perturbed values for the population numbers. Then, 
becomes just a given small number (which value we do not need to know here). Therefore, although the analytical 
form of the GFs, F 2 % IA ,Ffl IA , coincides with the F 22 , F^ , but the population numbers in the numerators, which 
determine the spectral weights, are slightly redistributed. Since a further analysis is possible only numerically, in 
order to get some impression of what is changed let us formally write down the population numbers in terms of O 
using the sum rule for the population numbers and the GFs gI 011 ' 111 !! 111 . 011 !, g[ioi,iii][ih,ioi] ) g [iio,iii][iii,ho]. 

^111 = i — TTTr i T> Who = 



l + 6[ei + e 2 ]' l + e[ei+e 2 ]' 

Qei 9e 2 
iV on — -i , ^ r ; t> ^ v ioi 



l + 6[ei + e 2 ]' l + 6[ei+e 2 ] 

ei = e^ u / 2 -»- h °ff\ e 2 = e ^ u / 2 -^. (B45) 

Since at X k — we have that P 2 = P 2 = iV° 10 + = N® 10 and equal to unity when the degeneracy is removed, it 
is clear that (1 — P 2 ) <C 1. Using this, we can write down an approximate expression for the Fermion GF, 

pHIA ^ 1 ~ p 2 , p 2 fB46) 

33 ^ + f + (l-P 2 )A fc w -f + A fe P 2 ' V ' 

We can neglect a weak and non-important dispersion oc (1 — P 2 )\k in the lower pole since this Hubbard sub band 
is fully filled and does not contribute to the cohesive energy of the system. However, the important thing following 
from this expression is that due to the derealization of the orbital in the region of energies uj ~ ^ the bandwidth of 
the band decreases, the spectral weight in upper pole is also decreasing from unity to P 2 < 1, and in a lower pole a 
small, but non-zero weight arises. 

Thus, we see that the function Fyy is reduced indeed to the form of F exp - t yy only in the case when P 2 = Nhq + 
Nm — 1. Since all other population numbers are not equal to zero, P 2 < 1 and the picture based on an effective 
single-electron potential should be corrected. 

Switching on mixing does not change the picture in this approximation, but slightly Shifts the value of <d and P 2 . 

What physical consequence can be expected of those shifts of the spectral weights? 

First, since a non-zero spectral weight of the orbitals 1 and 2 has appeared at the energy iuj ~ U /2, and the energetic 
barrier between the neighboring ions is not that wide nothing prevents a hopping between these orbitals belonging to 
different ions to develop, as well as mixing of these orbitals with the conduction bands. Therefore, the self-consistency 
loops should lead to the formation of bands due to hopping, etc. It should be noted, that the magnitude of the 
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shifted spectral weight is small because its value is determined by the competition between the Hund-rule intra- 
atomic exchange energy and cohesive energy due to a derealization of the empty part of the-;f-orbftals. The latter 
essentially depends on the degree of proximity of the upper transition A2 to the Fermi energycil, and, as will be seen 
below, the Hund-rule integrals split the upper transitions. This also leads to a suppression of the band formation 
with the transitions which are much above Fermi energy. 

Second, the expectation that pure Fermion bands will be developed from the "upper" orbitals happen to be valid 
within the approximation Nm — 0. Actually, there is no purely single-electron bands. 

Third, these bands are narrowed by the spectral weights. Although this is a well known effect from the studies 
of the Hubbard models, this is important in our discussion of an possible improvement of the ab mifia-calculations: 
even when SIC mimics the orbital polarization and improves essentially the description of the matteiO, it does not 
provide a fully correct physical picture. 

Fourth, not only the bandwidth of the band t vv (k) is decreased, t vv (k) t V y{k)[N lw + Nm], but, also, a part 
of the spectral weight is moved from the high-energy region iu> ~ to the low-energy iuj ~ Aj (see the Eqs for 
the Fermionic GFs F). Since usually the main contribution to the value of the moment comes from the localized 
orbitals, one can expect a slight deviation (decrease) of this moment from an integer value. In rare earths this effect 
is weak, but it is clear that, for example, in compounds with intermediate valence the deviation of the moment from 
integer is quite strong. Thus, this effect of the transfer of the spectral weight occurs due to derealization of one of the 
transitions (detailed consideration, of the mechanism of derealization due to a correlation-caused transfer of spectral 
weight is given in a recent paperEiJ). 

Fifth, when a long-range magnetic order arises in the system, the effective field, shifting the lower levels, as discussed 
above, changes the magnitudes of all population numbers (again due to the sum rule); the latter may change the 
bandwidth. Thus, if the correlated bands are involved into the formation of magnetism, the essential physics of the 
formation of the magnetic moment consists not in a rigid (self-consistent) shift of the bands, but in the changes of 
the bandwidths involved. From the technical point of view the less on w hich o ne ca n extract from the considerati on o f 
this model is the following. The speculations, based on Equations ( B22| ) and (B24) and which lead us to the Eq.(B25) 
might be considered as a ground for the ab initio calculations either within the standard model for the lanthanidcs 
and, also, within the SIC method, where SIC is applied to some of the orbitals. Some of the orbitals are fully localized 
and are placed very deep in energy either "by hands" in first case, or by a SIC potential i n the second one, while 
only the orbitals remaining at higher energies are delocalized. However, as we see from Eq. ( p346 ) this picture is not 
fully correct it misses the shifts of the spectral weights. As we shall see below, a very simple recipe to correct for this 
exists. In spite of the simplicity of the approximation, Hubbard-I, used in this example (which is the same in spirit as 
the ARF, discussed in the previous section), should reflect the mechanisms of the spectral weight transfer correctly 
since the driving mechanism is the sum rules, which comes from the exact commutation relations. Therefore, the 
mechanisms discussed above are expected to remain valid in higher approximations. 
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FIG. 1. The first contribution to the GF comes from the hopping/mixing matrix element. The solid line is the locator, T>, 
the circle is the spectral weights and the wavy line is the hopping/mixing. All lines in further graphs can be renormalized by 
this insertion. 



FIG. 2. The second diagram, a, appearing in SCPT, The dashed line is the Coulomb interaction. 
FIG. 3. The exchange graph, b, corresponding to the expression given in Eq. (|S6"|) . 

FIG. 4. A loop diagram in SCPT, graph c. The zig-zagged line represents the end in e, the constants of the algebra. The 
thick line represents the incoming and outgoing indices in e, i.e. 1" and 2 in e\ b 2 ■ The corresponding analytical expression is 
given by Eq.@. 

FIG. 5. Contribution d in Eq.(^). The curly line represent the Bosonic Green function (TZ 6b Z ib ). 

FIG. 6. Contribution e in Eq.(||). 
FIG. 7. Contribution f in Eq.@. 
FIG. 8. Graph which has an exact analogy in WCPT. A two-exchange process. Contribution g. 
FIG. 9. Graphs which has an exact analogy in WCPT. Contribution h. 
FIG. 10. Graphs which has an exact analogy in WCPT. Contribution k. 
FIG. 11. Graph which has an exact analogy in WCPT. An exchange diagram with a loop insertion. Contribution 1. 
FIG. 12. Graph which has an exact analogy in WCPT. An exchange diagram with two loop insertion. Contribution m. 

FIG. 13. Graph with one line renormalized by the first correction, Fig. |^. Any locator line can be dressed by the graph in 
Fig. ^. Contribution n. 

FIG. 14. The screening process, the localized transitions are screened by the delocalized transitions. The screening of the 
delocalized transitions by localized is very weak due to the large energy gap involved. 



36 



FIG. 15. The screening process can be continued, however, only delocalized transitions screen the Coulomb interaction. 
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